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EXECUTIVE SUMMARY 


Between 1984 and 1992 four major test series were performed in the HDR containment 
encompassing various fuels and three different axial positions in the high-rise, multi-level, 
multi-compartment facility. At that time, each HDR fire test series was accompanied by 
extensive efforts to evaluate the predictive capabilities of a variety of fire models and codes 
developed in different countries by both blind pre-test and open post-test computations. A quite 
large number of open issues remained in the area of fire computer code predictive qualities upon 
completion of the HDR program. 


In the meantime, large progress has been made in improving and consolidating fire models and 
computer codes of all levels of simulations. This progress merits revisiting both experimental 
results and fire computer code validations. The results of the research efforts for this grant 
during FY 1997/98 are documented in this volume: 


Volume 3: Test Series Description and CFAST Validation for HDR T51 Wood Crib Fire Test 
Series 


Volume 3 by focusing on the HDR T51 wood crib fire experiments covers the following aspects: 


e Section | provides an overall introduction to the HDR test facility and especially the 
containment building layout. It provides an overview of all four major HDR fire test groups 
utilizing a range of fire sources including: propane gas burners, wood cribs, liquid fuel pools 
and nozzle releases, and prototypical electrical cables. These fires have been set at three 
different axial elevations within the containment building under natural, forced, and 
combined ventilation conditions. 

e Section 2 gives a detailed account for the compartment layouts for the wood crib 
experiments. It also lists all fuel and thermophysical material properties involved in the 
experimental setup. 

¢ Section 3 describes the objectives, requirements, and functional principles of the 
instrumentation applied during the test series and documents the positions of all sensors used 
in both tabular and graphical forms. 

e Section 4 briefly summarizes the common test procedure used for executing every 
experiment. 

e Section 5 provides an overview of major experimental results of the wood crib tests in two 
subsections. First, selected transient histories are shown for temperatures, gas 
concentrations, and velocities in the different connected compartments, including the dome, 
for the three experiments spanning the range of gas fire powers examined. The second set of 
experimental results involves the maximum values of the same quantities as a function of the 
applied fire power. 

* Section 6 addresses numerous aspects of potential contributions of the wood crib 
experiments towards the validation of zone model codes such as CFAST (see Volume 2) and 
containment system codes such as GOTHIC. 

e Section 7 describes the CFAST model developed for the wood crib tests. 

* Section 8 discusses the results of computations using the CFAST model described in section 
7 for each of the three wood crib tests. 
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¢ Section 9 addresses the accomplishments achieved and deficiencies experienced with CFAST 
while modeling the wood crib tests. 

e Appendix A contains the CFAST input files for the model discussed in Section 7. 

¢ Appendix B contains a discussion of the spurious results that are obtained when trying to 
implement the CFAST O2 card while defining the fuel in an input file.. 
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1 INTRODUCTION 
1.1. HDR Test Facility and Containment Building 


The HDR (Heiss-Dampf Reaktor) facility, shown in Figure 1.1, was the containment building 
for a decommissioned, experimental reactor in Germany. The building, while smaller in volume 
than a typical US containment building, contained many features which made it valuable for use 
in a containment research program. Many of these features also make it extremely valuable as a 
generic source of test data for industrial facilities. The building was a cylinder approximately 20 
m in diameter by 50 m in height topped by a 10 m radius hemispherical dome for a total facility 
height of 60 m.. Internally the building was divided into eight levels with each level further 
subdivided into smaller compartments. For a typical HDR test approximately 60-70 
compartments were available. Compartments were connected by a variety of flow paths which 
included doorways, pipe runs, cable trays, hatches, and staircases. Three fixed and two 
adjustable vertical channels were provided for in the form of an elevator shaft, two staircases, 
and two sets of equipment hatches running the axial length of the building which could be 
opened or closed to change the available vertical flow path at each level. Much of the original 
equipment from the nuclear steam supply system was still present in the facility including the 
reactor vessel, primary and secondary piping, pumps, electrical connections, and ventilation and 
exhaust systems. The total free volume of the facility was 11,000 m* of which the dome 
contained 4,800 m°* above the operating deck. The HDR containment, its compartments, and 
internal structural materials, vent flow openings and other pertinent data are documented in [1]. 


Steel shell: 
Cyl. part: 30mm 
Hight: 60m 
Diam.: 20m 


Fire level 


Oil T52 amp Containment: 


Free volume: 11384m° 


L1.900 Compartments: 62 
L1.800 Intervention Concrete: 7700m° 
L1.700 area 

L1.600 


Fire level 
<jeom= E41.1-10; 42 Oil and cable 


<j 751.1-2 Gas and wood cribs 


Figure 1.1: HDR Facility and Fire Test Group Locations 
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1.2. Summary of Fire Test Matrix 


From 1984 to 1991 a total of four test series divided into seven fire test groups were performed 
inside the HDR facility. The fire tests consisted of the T51 series, six propane gas tests, three 
wood crib tests, and five more propane gas tests; the T52 series, four hydrocarbon oil pool tests; 
the E41 series, ten hydrocarbon oil pool tests; and the E42 series, three cable fire tests. Figure 
1.2 shows the overall test matrix and range of fires powers tested and Figure 1.1 shows the 
location of the various test series inside the HDR facility. Each test series was performed at a 
different location inside the containment building as indicated. 


10000 


kW 


5000 


1@) 4 
Gas 2 Hydrocarbon | Hydrocarbon | Hydrocarbon} Cable Fire 
Open Open Closed 
T51.11-19 ‘84/85 751.21-25 '86 |T52.11-15 '87 | E41.1-4 '88 E41.5-10 '90 E42 '91 


Level 1.400 Level 1.900 Level 1.500 


Figure 1.2: Fire Test Group Summary 
The fire tests were performed with the following general objectives: 


e An improvement in the general understanding of fire phenomena including smoke and 
aerosol production, distribution, and removal; temperature and pressure changes; and 
transient combustion in a large scale building. 

¢ A better understanding of the effects of boundary conditions on fire phenomena. 

¢ The creation of a large database for fire model and fire computer code validation. 

¢ An increase in the ability to plan for successful fire fighting and rescue operations inside a 
burning high-rise structure. 


The multi-level, multi-compartment structure of the HDR facility with its vertical shafts, large 
dome, and concrete and steel construction means that subsets of the fire test database have 
applications outside the nuclear industry. In general the fire test data can be used to gain insight 
on many industrial and commercial facilities as most share basic HDR features such as being a 
multilevel, steel and concrete structure with ventilation systems. More specifically, data from 
the large dome can be applied to hangars and atrium spaces. Data from the vertical shafts can be 
applied to any facility containing elevators, large vertical pipe channels, etc. 


Introduction 1-2 


T51 WOOD CRIB FIRE REPORT NUMAFIRE:03-98 


Each individual test series had its own specific objectives, which have been specified in the 
respective test series Design Report containing all pertinent geometric data, initial and boundary 
conditions, instrumentation plan, test procedures, and summary descriptions of the computer 
codes that participated in the pre-test and post-test computations. Data Reports were issued right 
after the experiments were performed and contained corrections/modifications of test 
procedures, qualification of the sensor operability and quality as well as all measured data in 
plots. All documented data have been stored on the PHDR data bank with the same format and 
sensor descriptions as used in all other HDR safety research experiments. Quick Look reports 
present and interpret the data according to the test series objectives and the associated physical 
phenomena. In addition to the presentation of the data of the individual experiments, results 
across the test series are documented. Moreover, Quick Look reports contain the comparisons 
between data and blind pre-test computational results by different models and codes used by the 
respective group of national and international participants. The Final Evaluation report 
documents all data assessments from the test series together with final conclusions and open 
issues. In addition, it contains the comparisons between data and open post-test predictions and 
identifies the learning effect, model and code improvements observed, lists remaining 
discrepancies, and open modeling issues. It is the final document for the test series. Section 7 
lists all relevant documentation cited above for the respective HDR fire test series. The 
respective reports will be referenced where applicable in Section 1.3, which summarizes the fire 
tests. 


The T51 test series, performed at the 1.400 level in the lower portion of the containment, was 
designed to be a relatively low power, exploratory test series in order to determine basic 
parameters of fire phenomena inside the facility [2-10]. The temperature changes inside of the 
fire room and the spread of smoke through the building and building ventilation systems was 
examined to determine safety margins for future, higher powered tests. 


The T52 test series, performed just below the operating deck, was designed to simulate a large 
cable fire through an equivalent oil fire [11-13]. The effects of ventilation systems on smoke 
movement was examined to assess rescue and fire fighting techniques. One major objective was 
to measure the plume behavior from the fire into the dome. 


The E41 test series, performed in the level above the one for the T51 test series, incorporated 
experiments that spanned the total range of fire powers examined in the HDR facility [14-20]. 
Additional parameters examined during this fire test series were the effects of opening and 
closing doors to the fire room, filter loading rates, and the effects of fire suppression systems. 


The final test series, E42, was performed at the same level as the E41 tests. The tests, consisting 
of cable fires, were to collect data on the burning of prototypical cables in cable trays under 
natural conditions [21-25]. The fires took place in an completely isolated set of 
subcompartments to prevent the spread of toxic combustion products, namely dioxin, resulting 
from the burning of the PVC insulation. A primary objective of these tests was to monitor the 
propagation of the fire through racks of cable trays in various orientations and to closely 
examine the spread and impact of combustion products. 
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Initially, the HDR fire tests were designed, performed, and evaluated solely by the Nuclear 
Center Kalrsruhe, German universities, industry, and research labs. However, the international 
nuclear community quickly realized the value of these tests [10]. Which resulted in international 
support, cooperation, and participation throughout much of the fire testing program at the HDR. 
Reflecting this is the fact that one of the E42 tests was selected to be a European Commission 
Standard Problem for the evaluation of computer fire models [24,25]. 


1.3 Overview of Individual Fire Test Series 


With the large variety of fire experiments performed in the HDR over many years, it is 
important to see where any one particular set of tests fits into the overall database of 
information. To this end a brief description of each of the fire test groups follows. 


1.3.1 @Gas\FireTPests|\@P51.UIFT SIC IS eis 19 ands P5122 151225) 


The gas fire tests, the T51 test series [2-10], were the first set of fire experiments performed in 
the HDR facility, and they are the subject of the remainder of this volume as well as volume 2 
[26]. A total of 14 tests were executed between 1984 and 1985. These tests consisted of three 
subgroups of five gas fires, a single gas fire performed at the end of the wood crib test series [5], 
and five additional gas fires [6-8,10]. The tests all took place in a specially constructed fire 
room on the 1.400 level, shown in Figure 1.3, of the HDR facility. This fire room was 
connected to a hallway which terminated under a vertical shaft formed by open maintenance 
hatches. Each experiment followed a similar test plan of a short period of pre-fire data 
collection to record initial conditions, followed by an hour long fire, and ending with 
approximately half an hour of cool down data collection. The fuel for each of the test was 
propane gas intended to be premixed with 10% excess air drawn from a vent in room 1.603. For 
the first group of gas tests no ventilation systems other than the air supply for the gas burners 
was employed. For the second group of gas tests a vent was constructed which connected the 
fire room to the 1.600 level. The vent had an adjustable damper which could be controlled 
during an experiment to change the size of the vent opening. 


This first test series had a number of primary objectives. The foremost objective was to 
demonstrate that fire tests could be performed safely inside the HDR containment building as the 
integrity of the structure was still regulated as a nuclear facility. Another objective was to 
determine the extent to which the fire would involve the building in its entirety. A further 
objective was to examine the ability of the ventilation systems to remove smoke and other fire 
products. Lastly, data collected during the tests would serve as a initial data for computer code 
evaluation. 
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The gas fire tests contain a number of characteristics which pose different of challenges for fire 
code models. These are: 


¢ The fire room is not a rectangular parallelepiped. The floor cross-section is L-shaped as can 
be seen in Figure 2.1 of Section 2. This geometric irregularity acts to impede some of the 
mixing that would otherwise occur in a symmetric compartment. 

¢ The fire source is not a single location on the floor in the center of the room. Rather, there 
are Six gas burners mounted on the wall 0.375 m off the floor along the L-side of the rooms 
length. Therefor the fire cannot be truly considered a point or local area source for the 
purpose of evaluating mixing and entrainment using common zone model approaches. Also 
the presence of the wall that the burners are mounted on prevents the formation of a typical, 
axi-symmetric plume that is assumed in many fire models. 

¢ The number and selection of burners used varied depending on fire power. 

¢ The doorway of the fire room is located at a corner, rather than at the center of one of the 
room’s walls. As with the shape of the room this affects the mixing that takes place inside 
the fire compartment. 

¢ The hallway from the fire room terminates in a subcompartment with a narrow vent, 0.5 m 
high, along the floor and a ceiling vent to a shaft leading to the hemispherical dome. 
Therefore, a fire model must be capable of handling a large ground level airflow as well as a 
separate, large buoyant plume in the same compartment. 

¢ The hallway from the fire room is not a rectangular parallelepiped. It is a volume of 
revolution, a rectangle slowly increasing in width rotated at a fixed distance about an axis. 


Table 1.1 on the next page contains a brief summary of the major characteristics of the gas fire 
tests. Figure 1.3 shows a top view of the fire floor. 
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Figure 1.3: Level 1.400, Fire Floor for the T51 Tests 
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Table 1.1: Gas Fire Test Series Summary 


Gas Ventilation and Other Test Burners 
Consumption Execution Comments Used 
(m*) 
TSN All run with the same 
Pit 12 14.63 configuration with only fire 
T51.13 26.62 power changing Burners 2,3,4,5 


between 1.600 and 1.700 
Uses Wood Crib sensor map 

dis iyd! 716 Pail wis Changes in sensor map Burners 1,2,5,6 
Repeat of test T51.13 with vent 
to 1.600 closed 

Aote2? eS 27S 30 minutes with vent 100% Burners 2,3,4,5 
open 
15 minutes with vent 75% open 
15 minutes with vent 25% open 
to 1.600 closed 

TS L24 951 36.58 30 minutes with vent 100% Burners 1,2,3,4,5 
open 
15 minutes with vent 75% open 
15 minutes with vent 25% open 

30 minutes with vent 100% Burners 1,2,3,4,5 


open 
1.3.2 Wood Crib Fire Tests (T51.16-T51.18) 


30 minutes with vent closed 


The wood crib tests were part of the T51 series of experiments [5,8,9]. The wood crib tests, 
while not a fuel typically available in a nuclear power plant, were added for the benefit of the 
fire community which does use wood cribs as a standard fire load. Three separate tests of 
increasing fire power were executed. The tests took place in the same fire room as the gas fire 
tests. Each test consisted of burning one or more cribs made up of 30 cm x 4 cm x 4 cm beams 
of pine containing 8% humidity. The beams were nailed together into 15 layers of 4 beams each 
with adjacent layers having a 90° rotation of the beams, Figure 1.4 shows the construction of a 
wood crib. A 300 ml reservoir of mineral spirits was used to start the ignition of the wood cribs 
which were allowed to burn uncontrolled. Electronic scales underneath the wood cribs were 
used to determine the time-dependent burning rate for use as input functions for the computer 
code simulations. As compared to propane gas which burns relatively smokeless, these wood 
crib tests were performed with a main purpose of evaluating the response of the HDR facility 
and ventilation systems to heavy loadings of smoke in an effort to determine safety margins for 
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future oil fires. The wood crib fires lasted on the order of 30 minutes. Table 1.2 gives some 
additional details on the wood crib tests. 


The wood crib tests produced large quantities of smoke which were quickly distributed 
throughout the whole containment. This smoke overloaded the building ventilation system’s 
HEPA filters and resulted in adjusting the testing schedule to accommodate the longer time 
required to clean the containment atmosphere between tests. The smoke was corrosive to the 
test equipment of other experiments, and some instrumentation was damaged. The smoke 
deposits of the HDR surfaces also proved difficult to remove, with success only occurring in 
cleaning of metal surfaces. 


Figure 1.4: Wood Crib Construction 


Table 1.2: Wood Crib Fire Test Series Summary 


Wood Ventilation and Other Test 
Consumption Execution Comments 


Start of Wood Crib sensor map. 
Fires were naturally ventilated 
and natural convection 
conditions existed in the 
containment. 


(7 cribs) 
(11 cribs) 


Introduction 1-7 


T51 WOOD CRIB FIRE REPORT NUMAFIRE:03-98 


1.3.3. Oil Fire Test Summary (T52) 


The second test series of fire experiments was the T52 oil fire test series which consisted of four 
oil pool fire tests performed in 1986 [10-14]. The tests ranged in power from two to four 
megawatts with the fire lasting approximately 30 minutes. Whereas the previous test series, the 
gas and wood fires, were performed at a level low in the containment building it was decided to 
position this test series high up in the containment building as shown in Figure 1.1. Thus, the 
fires were positioned in a special fire compartment constructed on the 1.900 level, the level just 
below the operating deck. It was anticipated that this would confine smoke and soot to the dome 
region. The fire compartment, shown in Figure 1.5, was located such that it vented directly into 
the dome through the maintenance hatch next to the spiral staircase. Fuel for the fires consisted 
of an initial volume of oil in a pool with a surface area ranging from 1 m’ to 3 m” in size. 


Figure 1.5: T52 Oil Fire Compartment 


The initial amount of fuel was augmented by a nozzle feeding a continuous supply of oil once 
the initial pool was consumed. Each fire lasted approximately 30 minutes. Oxygen for the fires 
was supplied either by natural convection alone or a combination of forced and natural 
convection. 


For this test series special attention was paid to the buoyant fire plume entering the upper dome. 
Two-dimensional grids of thermocouples and other sensors were placed at two axial levels 
within the plume to aid in determining the plume’s evolution in the dome. 


In addition to the generic purposes of improvements in knowledge about fire dynamics in a 
complex structure this test series introduced the concept of selective pressurization of test 
compartments for the prevention of smoke entry in rescue/intervetion areas. For this test series 
the elevator shaft next to the main staircase, see Figure 1.1, was pressurized and monitored to 
determine if selective pressurization was indeed capable of maintaining the entire shaft as a 


relatively smoke free area for the purpose of evacuation or for the staging of emergency 
personnel. 
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Some of the significant results are noted below: 


e The fires quickly reached flashover conditions, turning the fire room into a large fire ball 
with heavy soot production. 

¢ As the fire vented directly into the upper dome a large buoyant plume formed whose basic 
characteristics were measured. 

¢ The large buoyancy forces of the plume rising through the maintenance hatch behaved like a 
jet pump; that is large quantities of air were entrained into the plume which resulted in a 
large global circulation inside of the entire facility which widely spread the soot throughout 
the whole building. 

e Provided a sufficient air flow rate was used, the selective pressurization strategy was 
successful in keeping the elevator shaft free of smoke. 

¢ Due to the high entrainment, fire plume temperatures impinging on the containment steel 
shell were rather low. 


Table 1.3 below summarizes some details on the T52 tests. 


Table 1.3: T52 Oil Fire Test Series Summary 


Test Peak Fire | Pool Size | Initial Fuel | Fuel Delivery 
Power (m”) Volume Rate 
(kW) (liters) (liter/min) 


1.3.4 Oil Fire Test Summary (E41) 


The T52 test group indicated that both higher power and longer duration tests could be withstood 
by the HDR facility. A further set of oil fires, the E41 test group [14-20], was performed to take 
advantage of this. This test group, which consisted of ten tests ranging in power from six to ten 
megawatts, took place on the 1.500 level of the containment building. As with the other test 
groups a specially prepared fire compartment was used for this series. This compartment, shown 
in Figure 1.6, was significantly larger than compartments for the other tests and included 
sprinkler systems, ventilation systems, and a remotely operated doorway. For this test series the 
building ventilation systems were equipped with different types of filter setups. Furthermore, 
autonomous, aerosol measurement devices were added to the sensor equipment. 


The addition of extra features to the fire room and ventilation system allowed the examination of 
some additional fire phenomena. Filter loading and clogging was examined through the use of 
the different filter systems. The effects of steam release into the fire room was examined. The 
interrelationships of doorway openings and mechanical ventilation were explored. The selective 
pressurization strategy was examined further. Tables 1.4 and 1.5 provide details on this test 
group. Note that each test in the latter portion of this test series actually consists of a series of 
individual subtests. 


Introduction 1-9 


T51 WOOD CRIB FIRE REPORT NUMAFIRE:03-98 


YN 
NK ee 
SNS Ni 
tiem [RSI \A\ 
ANY 
Curtain HHH AVX 


SA > 
1.504 SN 
+2.8m 


Figure 1.6: E41 Oil Fire Compartment 
Some of the significant results of this test group are given below: 


e Fire extinguishing systems were tested under extreme conditions of fire power and 
temperature due to the high fire powers, as high as 10 MW, in the fire compartment. 

¢ Spatial and temporal distributions of aerosols were measured at different locations. 

¢ Depending on the ventilation system settings a variety of flow circulation modes were 
observed inside the containment building. 

¢ Selective pressurization of the elevator shaft was again successful in preventing smoke from 
entering this rescue shaft. 

e Filters continued to become overloaded with soot even when a prefiltered bank consisting of 
coarse filters was added to the filtration system. 


Table 1.4: E41.1-10 Oil Fire Test Series Summary 


Pool Size (m7) and Max Power 
Pool Wall Material (1) (kW) (min) 


"For tests E41.5-10 the fuel volume represents the initial pool volume. 


E41.5 2 (steel) 
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Table 1.5: E41.5-10 Oil Fire Test Subsection Summary 


En a 
Test Subsection Time Fuel Addition Door 
(kg/s) 
(min) 


F41:5ia 5 Initial Volume | Closed 
E41.51b 
BEALS E41.52 
1 


aaa sirtt a 
8a 
Fan hd 


d 
a) 
2) 
5 
2) 
5 


oe Door 1.45" 
E41.64 


cone 
PP E4166 [80 


E41.7 


Be 6 90 ook |e .03 4055... Doorie4s*s| 

Both Open 
res 
| £41.84 | 60 |  0.03-0.05 | Door 1 Open | 
PRATER Oman emeOOM MIMD TP None? Siclosed 0 F" 
~ 
vee ee 14 Ome | 05-0.079 | Closed ener 
EA cee to ae None 028 | Closed 
Initial Volume 
0.05 
0.03 
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1.3.5 Cable Fire Test Summary (E42) 


The cable fire test group was the last set of fire experiments performed in the HDR and had the 
primary purpose of evaluating the effects of a prototypical fire using real fuel sources, e.g. the 
electric power and instrumentation cables used in power plants [21-25]. Due to concerns of 
dioxin production from the PVC cable insulation, this test group was performed in an isolated 
subset of compartments on the 1.500 level which is shown in Figure 1.7. Additional partitions 
and ventilation and fire extinguishing systems were constructed on this level to prevent the 
spread of toxic combustion products through the rest of the facility and into the local 
environment. Three tests involving different amounts and types of cables were performed. It is 
important to note that the fire compartments were completely sealed for the duration of this test 
series which created problems in determining the exact fuel source available or consumed during 
any given test. As shown in Figure 1.8, before the first test, E42.1, many of the cable trays were 
wrapped in Alsiflex mats in an attempt to prevent the combustion of those cables during the first 
test. Attempts were made to isolate specific cable trays from burning by covering some of the 
cable trays in Alsiflex blankets which could be removed for other tests. The blankets did not 
completely prevent combustion of the protected cables; that plus a lack of information on the 
fraction of exposed cables which completely burned results in an uncertainty in specifying the 
exact fuel source available and consumed during each test. 
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Figure 1.7: E42 Cable Fire Room Figure 1.8: Cable Tray Layout 
Some of the important results from the E42 test series are given below: 


¢ Depending on the particular configuration of available cables the cables fires were either self 
sustaining to the point of flashover or burned out after a short period of time. 

¢ Dioxin production from the PVC insulation was not detectable/measureable. 

¢ The fires were capable of becoming intense enough to burn the cables underneath the 
Alsiflex blankets. 

e The presence of the blankets actually acted to prolong fires as they prevented water from the 
sprinklers from reaching the cables under the blankets. 
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2 FACILITY DESCRIPTION 
2.1 Compartment Layouts for the T51 Wood Fire Tests 
2.1.1 Fire Floor (Level 1.400) 


To avoid damaging the structure of the HDR facility, as the building was still considered a nuclear 
facility, a special set of fire test rooms was prepared at the 1.400 level of the containment 
building, see Figure 1.1 for the location of the 1.400 level. These rooms also served to control 
the flow of gases in and out of the fire room. Figure 2.1 shows a cross section view of the 1.400 
level and indicates the location of these rooms which consisted of the fire room with a narrow 
doorway, a long hallway wrapping around the reactor vessel shield wall, and a curtained area 
centered beneath the maintenance hatch next to the main staircase. For the remainder of the 
facility no special precautions were undertaken with respect to insulation as gas temperatures 
outside the fire floor were anticipated to be below damage causing levels. Table 2.1 below gives 
the geometric data of the prepared compartments [2]. These compartments are the same 
compartments that were used for the gas fire tests T51.11-15, T51.19, and T51.21-25. 


Table 2.1: Fire Compartment Dimensions 


Compartment | Height Area Volume Doorway 
(m) (m2) (m3) (m wide x m 


| FireRoom | 2.750 | 9.66 | 26.58 | 1.01x1.975_| 


N/A 


The fire room, its vertical cross section through the 0-180° plane and horizontal cross section at 
the +0.0 m elevation are shown in Figure 2.2, was constructed inside of room 1.405. The floor, 
walls, and ceiling of the fire room were lined with 25 cm of Ytong fire brick. The ceiling, which 
would be exposed directly to the fire plume, had additional protection in the form of a 3 cm thick 
layer of Alsiflex fireproof matting. Alsiflex is a ceramic fabric. Due to damage from the first 
series of gas fire tests, T51.11-15, this layer was replaced before the start of the wood crib tests. 
The floor of room 1.405 below the Ytong fire bricks consisted of a 1 m thick layer of concrete 
coated with a 1.5 mm thick coating of paint. The wood cribs were placed in the larger open area 
just inside the fire room doorway, see Figure 2.2. The wood cribs were placed on top of a 0.281 
m high platform in the manner indicated in Figure 2.2. Lastly, on the ceiling of the fire room, over 
the wood crib platform, a fire suppression steel pipe was installed. This pipe would inject coolant 
water to extinguish the wood crib fires in the event that unsafe conditions developed during a 
wood crib test. 
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Figure 2.1: Fire Floor at Level 1.400 for T51 Wood Crib Tests 
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Figure 2.2: Fire Room and Doorway 


The wood crib platform was designed to allow the measuring of the wood’s combustion rate by 
the use of three digital scales. The platform consisted of a scale plate with three layers of 
insulation to protect the digital scales from the high temperatures of the fire. Two layers of 
Promalan HT 400 with an intermediate layer of Promasil 1060 where used to protect the surface 
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of the scale plate, see Figure 2.3. The sides of the platform were protected by a wall of Ytong 
firebrick. On the fire room floor beneath the scale plate, two flat coolers (0.800 m x 0.826 m) 
were installed for further protection of the scales, see Figure 2.4. These coolers were cooled by 
chilled water. The whole platform rested on top of three digital scales which were used to 
determine the burning rate. 


Plate Cooler 


2.18m et! 


Figure 2.3: Wood Platform Cross-Sectional View Figure 2.4: Wood Platform Top View 


The doorway to the fire room had the same construction as the fire room itself. The doorway’s 
height was less than both that of the fire room and the hallway as depicted in Figure 2.2. 


The hallway, its top view shown in Figure 2.5, was constructed inside of room 1.406. This 
compartment wrapped around the reactor vessel shield wall, thus, the walls along the length are 
not parallel planes, but instead are concentric arcs. The doorway from the fire room enters the 
hallway slightly offset from the end of the hall. The walls and ceiling of the hallway were 
constructed of a 10 cm thick layer of Ytong fire brick. The floor of the hallway was not insulated 
and was of the same concrete and paint as in the fire room. This results in the floor of the hallway 
being 25 cm below the floor of the fire room and doorway. The hallway had a constant ceiling 
height, and its width varied from 1.14 m at the door to the fire room to 1.8 m at the entrance to 
the curtained area. This variation results from the curvature of the hallway and the general layout 
of neighboring compartments as depicted in Figure 2.5.. 


Hallway 


Figure 2.5: Fire Level Hallway Figure 2.6: Curtained Area 
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The curtained area, shown in Figure 2.6, was constructed next to the main staircase under an 
open maintenance hatch and was created primarily to force flue gasses from the fire away from 
the inside steel sheli and up through the vertical shaft of the open maintenance hatches by 
buoyancy. This area was created by draping Alsiflex mats from the ceiling down to 0.5 m above 
the floor. This 0.5 m gap served to enhance the cold air return flow through the hallway to the 
fire room and to force a separation of the hot and cold gas streams. It is important to note that 
the gap below the curtain opened up towards both the main staircase and towards compartment 
1.403 which had a floor elevation that was 1.1 m above the floor of the fire hallway or 0.6 m 
above the top edge of the gap beneath the curtain. The floor of the compartment below the hatch 
consisted of a 0.6 m thick layer of concrete protected by a 1.5 mm thick paint coating. The 
ceiling around the open maintenance hatch consisted of a 0.25 m thick layer of concrete protected 
with a 1.5 mm thick paint coating. 


Figure 2.7 shows a perspective view of all the added enclosures indicating (1) the fire room, (2) 


the doorway, (3) the hallway, (4) the area inside the curtain, (5), the curtain, (6) the maintenance 
hatch, and (8) a window for viewing the fire with a video camera. 
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Figure 2.7: Perspective View of Level 1.400 Fire Compartments 
2.1.2 Facility Remainder 
The following two tables, Tables 2.2 [1,26-27] and 2.3 [1,26-28], document the volumes of the 
different compartments in the HDR facility as well as the sizes of the major room interconnections 


available during the gas fire tests. Details of the layout of the HDR compartments can be located 
on the instrumentation maps shown in Section 3. 
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Table 2.2: HDR Compartment Volumes 


Compartment | Volume | Elevation | Height Comments 
Number (m*) (m) (m) 


8.50 


1.202 78 -9.20 5.80 Separate volume information not 
1.203 given for these compartments. 
1.303 


-5.80 


3 
-5.80 
-5.80 
a) Ws Separate volume information not 
ea | given for these compartments. 
Main staircase level 1.300 
Main staircase level 1.400 
Main staircase level 1.600 
Main staircase level 1.700 
Main staircase level 1.800 
Main staircase level 1.900 
| Sittiieunr Zhi iii PEs Not open during wood crib tests. 


1.402 40 3.50 No volume found, estimated from 
floor area. 


: 
Includes volume of fire room 


Includes volume of fire hallway and 
curtained area. 

Elevator shaft. Not open for T51 

158 4.50 4.50 | Separate volume information not 

pee be given for these compartments. 

eee eso ee ris ee 
eo: er ee ee CO eons Pert 
Oe ae eee 
aU Soe eee OMe, So if 
SE SRE ee eae ae 
CTS il ee Ee ee ee 
p53 8s | 3.50 | __ 8.80 _| Not open during wood crib tests. | 


Facility Description 2-5 


T51 WOOD CRIB FIRE REPORT NUMAFIRE:03-98 


a 

Number (m*) (m) (m) 
RBS Lee ie Ma 
ROR ee eS en 
[W604 se ist cea Se |NRAEOOO | MRS 25 ay: TSW npg ney a 
a a 
[1,606 7 2)" 8183 E0008 | ca 60 es a ee ne 


| 
i 
1 


Separate volume information not | 

1.608 siven for these compartments. | 
| 1609 | 59 | 10.00 | __4.75_| Not open during wood crib tests. _| 
192" |g 0,00: 094.75. |) he 
BO ee a a 
1.7010 8%. * oes Pe 016 Ot EE? 50h ae ee ee 
io de 203 T | IRS ATS OSES 2 OE ee ee 
805 14.25 Separate volume information not | 

1.901 siven for these compartments. | 


1.706 15.05 Not open during wood crib tests. 
9 
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= 
N 
Nn 
: 
— 
: 
= 
Nn 
: 
N 
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| __ 1, 707A) ee ATS IOS) 14.20 | ON i ai 
1. 70800 [it 590 MIANTSIOS IS (022535. OL 
. 1:80.10 fan) WS 43h TAO TOS. MRO. 80 5. IAM) Os 
|< WL 8O2HnhiNels f925 Wel HN20.60/ TH) GhTIOT ROO a) ON 
ja 1,804, | S798 E0160. ob TS00 | iis 
sl. 902.0 2s | 0190 825.30 AS OM he 
BE en en Oe ee ae eS 

1.904 Separate volume information not | 

1.905 given for these compartments. | 


1.906 25.30. (e450. Oils. | 0 ea 
94153 30.85 Cylindrical portion of dome 
2,660 40.00 10.00 | Hemispherical portion of dome | 
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2.2. Thermophysical Material Properties 
2.2.1 Thermophysical Wall Surfaces Properties 


There were five different materials which were used as compartment surfaces within the HDR 
facility. Alsiflex mats and Ytong firebrick were used to create the fire room and hallway on the 
1.400 level. In general, rooms in the HDR facility had painted concrete for the room surfaces 
with a different paint used for the floor than was used on the other room surfaces. Tables 2.4 and 
2.5 provides the thermophysical properties for these materials [16,20,35] 


Table 2.4 :Material Properties for Room Surfaces 


(kg/m’) (kJ/kg K) (W/m K) 


Table 2.5 : Thermal Conductivities for Room Surfaces 


Thermal Conductivity 
Material (W/m K) 


100°C _| 300°C 
Ytong Fire Brick | 0.09 | 0.15 [| 0.19 | 0.23 | 0.24 | 


2.2.2 Thermophysical Fuel Properties 


The fuel for this fire consisted of either five, seven, or eleven wood cribs constructed of pine with 
an initial humidity of 8%. Each crib had fifteen layers of four 30 cm x 4 cm x 4cm boards. Each 
layer was nailed crosswise to the layer below it, see Figure 2.8. The lowest layer of the wood 
cribs contained only the two side boards. The cribs were arranged on the raised platform as 
indicated in Figure 2.9. The crib designated number 2 was the ignition point, by means of a 
metallic can filled with 300 ml of alcohol, which was manually ignited for each of the three wood 
crib tests. Posttest data analysis used a heat of combustion value of 12600 kJ/kg to determine the 
maximum fire powers achieved [10]. 


Facility Description 2-11 


T51 WOOD CRIB FIRE TEST DESCRIPTION REPORT 


NUMAFIRE:03-98 


Figure 2.8: Wood Crib Construction 


Facility Description 


Figure 2.9: Wood Crib Placement 
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3 INSTRUMENTATION LAYOUT 
3.1 Introduction 


Because the fire research experiments were added to the HDR Safety Program about midway in 
the course of numerous safety experiments, the development of an instrumentation plan and the 
selection of the sensor types rested upon tested and proven measurement technologies. Most of 
these technologies were successfully applied during the previous HDR containment experiments. 
This proved to hold for the majority of typical pressure and temperature sensors. However, it 

_ was apparent from the outset that the fire experiments had somewhat different instrumentation 
criteria owing to the high temperature, low flow, and corrosive environment that the sensors 
would be exposed to. 


Therefore, the T51.1 experiments opened up new challenges for the instrumentation, especially 
because all previous expertise resulted from high-momentum driven flows only. On the other 
hand, the fire experiments required reliable instrumentation for buoyancy driven flows with 
much lower velocities and much higher temperatures. 


Given all these circumstances, test series T51.1 served as an exploratory test bed for advanced 
instrumentation such as velocity sensors, gas concentration sensors, smoke detection sensors, and 
other fire related sensors. 

In addition to the sensor qualification issues, numerous questions arose regarding: 


¢ safety procedures (injecting explosive gases in the building and the effects of high 
temperature loads on the structure) 
¢ the optimal placement of a limited number of sensors and sensor types. 


Answers and resolution guidance for both issues for all previous HDR experiments commonly 
rested on so-called design computations representing a broad spectrum of different approaches 
and models plus compliance with industrial codes and regulatory standards. 


In the case of the T51 fire experiments, all analytical and computational methods at that time 
were limited to treat only single compartment, single burning object and single vent flow 
opening. The resultant predictions were naturally overly conservative because they did not 
account for multi-compartment geometry, counter-current flows of hot and cold gas streams and 
associated mixing as well as heat transfer to structures. 


Therefore, safety measures and experimental procedures were extremely stringent and 
conservative, such as the installation of the curtain at the end of the hallway in order to keep the 
hot flue gas layer away from the inside surface of the containment steel shell. 


Equally, the instrumentation plan for T51.1 was primarily geared towards safety rather than 
towards physical phenomena in the first place. These details should be kept in mind when 
reading the subjects about sensor types described in Section 3.3 and the instrumentation map 
described in Tables 3.1 with its accompanying figures of facility cross-sections and instrument 
positions. 
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3.2. Objectives and Requirements 


With the background knowledge from the previous section, the instrumentation for the T51.1 
wood crib tests was designed to encompass the following elements: 


1. Instrumentation in the Fire Room: 
¢ 20 thermocouples in the fire room and door vent opening 
¢ 10 thermocouples at/in the structures 
* pressure sensor in the fire room 
e heat transfer instrumentation in the fire room and neighboring compartments 
¢ digital scales for measuring the combustion rate 


2. Instrumentation in the Containment: 
e thermocouples and pressure sensor from previous containment blowdown experiments 
¢ several flue gas analysis sensors 
e heat transfer blocks from previous containment blowdown experiments 
* pitot tube sensors 
e determination of smoke density 
e miscellaneous special sensors as describe in Section 3.3 


3. Instrumentation of Exhaust 
¢ thermocouples and velocity sensors 


4. Safety Instrumentation for Protection of the Containment’s Steel Shell Integrity 
¢ thermocouples along the height of the vertical staircase/maintenance shafts 
¢ thermocouples in the reactor dome above the operating deck 
¢ thermocouples at inside and outside steel surfaces 
* gas concentration sensors in the dome and lower containment regions 


Whereas the instrumentation in the fire room needs to satisfy the special fire requirements, the 
rest of the containment instrumentation relied upon the available, proven containment 
measurement sensors. All data acquisition needs were accomplished by the central HDR 
computer and data acquisition storage system. 


3.3. Instrumentation Descriptions 
3.3.1 Temperature Measurement 


NiCr-Ni, sheathed thermocouples were used for temperature measurements in accordance with 
German DIN 43710. The thermocouple sheath had a 3 mm diameter and an insulated tip. The 
signal wires did not require special treatment as long as they remained outside of the hot flue 
gasses. Depending on the thermocouples physical location in the facility, the signal wires were 
up to 20 m in length. As high frequency temperature changes were not anticipated, 
thermocouples with standard response characteristics were chosen; e.g. for temperature an error 
of +1% of the measured value and for strong thermal radiation conditions an error of +5%. 
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3.3.2 Pressure Measurement 


Figure 3.1 shows schematically the major elements of determining the pressure difference with 
the TELEPERM measurement converter. The difference between containment inside and 
outside pressures acts on the bellow and is transmitted through a lever to the flexible beam tube 
which in turn transmits to a differential capacitor providing an analog signal. The TELEPERM 
converter works for a pressure difference of up to 5 mbar with a response time of 0.3 s and a 
measurement accuracy of +1%. This device had to be protected from high temperatures; hence, 
its placement on Level 1.6 of the facility. 


E | 
SSIS Sa PSEN NAN SL 
Ne ae NS i 
#40 


Cable Protection 
Material 


10 :) 8 if Test 


2 


1 
1. Differential Capacitor 5. Bellow ’ 9. Phase Dependent Converter 
2. Control Lever 6. 10kHz Oscillator 10. Amplifier 
3. Flexible Connector 7. Voltage/Current UH Aux. Power Sensor Material 
4. Displacement Limiter Converter Ub_ Internal Power 
1A Output Signal 


Figure 3.1: TELEPERM Transmitter Figure 3.2: Local Heat Transfer Measurement 
3.3.3 Heat Transfer Measurement 


The local heat transfer at the containment steel shell was determined using the sensor depicted in 
Figure 3.2. The sensor used a 40 mm diameter annular control volume with a thin disk bottom 
of known material properties. Two NiCr-Ni sheathed thermocouples with a 0.5 diameter 
measured the disk temperature. 


Additionally, large concrete blocks were devised and equipped with thermocouples as 
schematically shown in Figure 3.3. The type of concrete chosen was the same as used for the 
construction of the HDR containment Except for the front surface of the block all other surfaces 
were insulated. These massive concrete blocks were positioned at location where high 
convective flows, such as in the staircases, could be anticipated. 


As the determination of the heat flux and subsequently the heat transfer coefficient at the 
measurement block’s front surface rests on the solution of the inverse heat conduction problem, 
errors in these quantities became larger when temperature differences between thermocouples 
became smaller. Therefore, the expected accuracy of these blocks was only +20%. 
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Figure 3.3: Heat Transfer Block 
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Figure 3.4: Gas Volume Analyzer 


One of the major overall objectives of the HDR fire test series was the evaluation of the hazard 
potential to personnel, fire fighting, and rescue teams dependent on the type of burning 
substance, ventilating conditions, and fire location within the high-rise, containment building. 
Aside from direct exposure to heat, it is the smoke and flue gas mixture (O,, CO, CO», C,Hm, 
NO,, SOx) as well as the production of HCI and potentially dioxin in the case of burning PVC 
cables which determines the hazard level. Therefore, instrumentation measuring the 
concentrations of these individual components had to be in place. Fortunately, as the majority of 
the experiments of the T51 test series used non-sooting gas flames, the requirements for smoke 
and gas analysis largely reduced to measuring O2, CO, and CO> concentrations. 


Figure 3.4 shows a schematic of the non-dispersive infrared photometer which worked with a 
modulated, single beam. This instrument allowed for continuos operation using a suction pump 
in the range of 10-100 1I/h volumetric flow. The device outputs a 0-10 VDC signal proportional 
to the volumetric concentration in terms of vol. % or ppm. 


Prior to the start of each experiment, these sensors were calibrated with a calibration gas. The 
measurement accuracy of these sensors was expected to be +2%. 


3.3.5 Optical Smoke Density (Extinction Coefficient) 


In order to follow the distribution and propagation of the flue gases inside the containment, an 
optical smoke densitometer, type ME82 made by Maurer, was positioned throughout the 
containment. A schematic of this sensor is shown in Figure 3.5, This sensor was used to 
determine the optical gas density in the rescue paths as well as the smoke density according to 
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German Standard DIN 4102 Pt. 1. The output from this device was converted to an extinction 
coefficient prior to data recording. 


As shown in Figure 3.5, a standardized light source in accordance with DIN 5033 emits a beam 
of light which passes through a control volume containing the gas to be analyzed. The control 
volume size can be modified. The amount of light passing through the volume is converted to 
an analog signal from 0-10 VDC corresponding to 100-0% transmittance. The measured values 
had an accuracy of +2%. 


1. Light Source 4. Smoke Channel 

2. Lens 5. Receiver/Filter TELEPERM 

3. Window 6. Output Signal See Fig 3.1 
Figure 3.5: Smoke/Gas Density Sensor Figure 3.6: Pitot Tube Velocity Sensor 


3.3.6 Velocity Measurement 


In order to measure the flow velocities in different regions inside the containment, pitot tube 
type sensors, shown in Figure 3.6, were used. These sensors determine the pressure difference 
relative to stagnation pressure and use that to calculate the velocity. Hence, they used to same 
TELEPERM transmitter as discussed in Section 3.3.1, Figure 3.1. To obtain the velocity, the 
gas density must be measured simultaneously. For this purpose a thermocouple close to the pitot 
tube was used to determine density by application of the real gas law. The pitot tubes were 
capable of operating in temperatures as high as 800 °C. As before, the TELEPERM transmitter 
was shielded to protect it from high temperatures. 


3.3.7 Video System 

3.3.7.1 Introduction 

The HDR facility was equipped with a color video system consisting of cameras, monitors, and 
tape machines. This system was used for monitoring the fire behavior in the fire compartments. 
The camera position for the T51 test series is shown in Figure 2.5. 

In addition, a black and white video network consisting of 20 cameras with a switching board 


was installed. This system was developed by the Technical University of Karlsruhe, Germany, 
for use in monitoring the evacuation of personnel from high-rise buildings during fire exercises. 
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3.3.7.2 Black & White Video System Network 


Figure 3.7 shows the black and white video network which was used to monitor smoke 
movement at up to 20 locations under low lighting conditions, 20 lux. The cameras are 
connected with 50 m long cables to a video switching board. This device switches to next 
camera after three half pictures are taken. With a camera frequency of 50 frames per second, the 
switching board could rotate through the cameras in 1.2 s. This results in a nearly simultaneous 
observation of the smoke throughout the building. The other elements shown in the Figure are 
self-explanatory. During test T51.18, the underventilated nature of the fire resulted in such high 


levels of smoke, that the video system optics quickly became coated in soot, rendering the test 


footage for T51.18 useless. 
: Video Mixer 


Camera Switching Device Bee 


Clock with 
Characteristic Data 


Tape Machine 


Montior 


Microphone 


Figure 3.7: B&W Video System 


3.3.8 Digital Scales 


The wood crib platform was supported by three digital scales, see Figures 2.3 and 2.4 in Section 
2.1. The scales were tared with the weight of the platform. Each of the three scales measured 
the force in Newtons exterted by the wood on the platform. To generate the burning rate the 
total of the three scales was taken and converted to kilograms of mass. 


3.3.9 Safety Measures 


In addition to the thermal insulations listed in Section 3.1, the measures described below were 
implemented for safety purposes: 


¢ The color video system monitoring the fire compartment was continuously operated. 

¢ A number of gas detectors were positioned at Levels 1.4 and 1.6 for safety reasons. 

e All other containment regions including the steel shell were monitored with thermocouples. 
e A fire suppression system was installed along the ceiling of the fire compartment. 
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3.4 Instrumentation Layout for T51 Wood Crib Tests 


This subsection describes the instrumentation mapping for the T51 wood crib tests. A complete 
listing of all instruments as well as diagrams showing their locations within the facility are 
documented in the tables following. To aid in reading the tables and diagrams the following 
nomenclature, standard for all HDR tests, is used for the instrumentation: 


CA: Force Sensor 

Clee Velocity Sensor 

GG: Gas Concentration Sensor 

CE Pressure 

(GOR Heat Transfer Measurement Block Sensor 
CS: Temperature Sensor 

i Temperature Sensor 

CV Velocity Sensor 

OA: Steel Shell Expansion 


In addition to the directly measured quantities, post processing was performed for some of the 
tests to yield indirectly measured parameters such as density and mass flow rate. These indirect 
measurements were not performed consistently throughout the test series. These measurements 
used the following nomenclature: 


CD: Calculated Density 
CM: Calculated Mass Flow Rate 
CQ: Calculated Heat Flux or Heat Transfer Coefficient 


In the table that follows sensor location refers to one of two coordinate systems. For heat 
transfer measurement blocks the location uses the front, center of the measurement block for the 
reference location with the position given in Cartesian coordinates [9]. All other sensors use the 
HDR center line at the +0.0 m elevation, see Figure 1.1, for the reference location with the 
position given in cylindrical coordinates [1]. 
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Table 3.1 lists all sensors in place for tests T51.16-T51.18. The table shows the 
quantity/parameter measured, the sampling frequency, and the location for each sensor relative 
to the appropriate coordinate system. Any special comments about the sensor’s performance is 
also given. Figures 3.1 through 3.14 schematically depict the sensors’ locations level by level in 
the HDR facility for Table 3.1. 


Table 3.1: T51.16-18 Instrument Network 


| 
Sensor Parameter Pasoucel ] Z Comments 
(HZ) | (cm) | (deg) | (cm) 
Pia Ric ad 
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ven ort dil ao ioeh Sl 
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— | fe he) 
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Table 3.1: T51.16-18 Instrument Network (continued) 


(Hz) 
(cm) | (deg) | (cm) 
1180 
450 
JR Sei ee ae ea) 
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7c AC) TT 
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i PMRO | 18801 2200 
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Table 3.1: T51.16-18 Instrument Network (continued) 
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Table 3.1: T51.16-18 Instrument Network (continued) 


Location 
Sensor Parameter Comments 

(deg) | (cm) 
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Table 3.1: T51.16-18 Instrument Network (continued) 


: 
Sensor Parameter oe cy R ) z Comments 
(cm) | (deg) | (cm) 
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Figure 3.8: Level 1.200 and 1.300 at the -6.5 m Elevation 
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Figure 3.9: Fire Room at the -0.7 m Elevation 
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Figure 3.10: Fire Room at the +0.2 m Elevation 
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Figure 3.11: Fire Room at the +1.1 m Elevation 
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Figure 3.12: Fire Room at the +1.7 m Elevation 
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Figure 3.13: Fire Room Vertical Cross-Section 
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Figure 3.15: Level 1.400 Velocity Sensors and Gas Sensors 
at the +0.0 m Elevation 
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Figure 3.16: Level 1.500 at the +6.0 m Elevation 
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Figure 3.17: Level 1.600 at the +12.0 m Elevation 
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Figure 3.20: Dome Level at the +31.0 m Elevation 
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Figure 3.21: Dome Vertical Cross-Section 
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4 TEST EXECUTION 


Test execution procedures for the wood crib tests were slightly more complex than for the gas 
fire tests. This extra complexity results mainly from pretest procedures for the handling of the 
wood. The procedures can be divided into pretest actions, test actions, and post test actions. 


Before any given wood crib fire test was executed, the wood was pretreated to in order to obtain 
an initial 8% humidity. 


The test execution procedure for the wood crib tests was as follows. Approximately one half 
hour before the start of each test the data acquisition system was activated. After verifying the 
operation of the data acquisition system the wood cribs were ignited. This was done by igniting 
a 300 ml of alcohol in a small canister located under crib number 2, see Figure 2.2. The ignition 
was done manually. The remainder of the wood cribs ignited from the radiative heat resulting 
from the ongoing fire. The wood crib fires lasted approximately 45 minutes. Data collected 
during each test was archived on reel-to-reel magnetic tape. During the archiving process the 
start of the fire was set to be zero minutes and ten minutes of data were stored from the prefire 
portion of the collected data . 


Following the end of a wood crib test, the HDR containment building was ventilated to purge it 
of smoke and combustion product gasses. The ventilated gasses passed through filters, and data 
was collected on the filter loading. The levels of smoke and soot produced during the tests, 
especially during the T51.18 test, required the complete replacement of filters after each test in 
order to keep filter system pressure drops below regulatory requirements. After ventilation 
further sampling was performed on the containment building atmosphere. Levels of hazardous 
gasses were checked to determine if further ventilation was required before personnel could 
enter the building. 


Following each test, cleaning was performed of the various sensors to remove any soot buildup 
that could interfere with their future performance. Test T51.18 with its high levels of soot 
production required a significant expenditure of manpower and materials to clean the facility and 
its sensors. 
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5 OVERVIEW OF EXPERIMENTAL RESULTS 


This section contains selected results from the wood crib fire tests performed in the HDR 
facility. Data from selected instruments for the three tests are shown in the first subsection to 
give a general overview of the transient histories and their similarities and differences. These 
are followed by results representing zonal values from the test series plotted as a function of the 
number of wood cribs at 16 minutes from the start of the wood crib fires. Plot legends in the 
second section indicate the instrument names and elevations; their positions can be identified in 
Section 3 of this report. 


ra | Selected Results 
Results from each wood crib test are shown in this section, tests T51.16 through T51.18. These 


tests involved the combustion of 5, 7 and 11 wood cribs respectively. The first figure, Figure 
5.1, shows the pyrolysis rate for each of the three tests. 
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Figure 5.1: Wood Crib Pyrolysis Rate (CA4500) 
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The following observations hold: 


Each test has the same basic trends in pyrolysis rate of a slow initiation phase lasting for 
about 10 minutes which is followed by a steep increase. Two progressively smaller peaks 
follow after the maximum pyrolysis rate, which occurs at 16 minutes, although in tests 
T51.17 and T51.18 there is overlap between these secondary peaks. 

The duration of significant combustion activity is nearly equal for all three tests, 
approximately 35 minutes. 

There is an increase in the time to peak pyrolysis between tests T51.16 and T51.17, but not 
between tests T51.17 and T51.18, and all three achieve a maximum pyrolysis rate at 16 
minutes. 

The increase in the number of wood cribs results in the increase of the maximum pyrolysis 
rate. 

Towards the end of the wood crib fire, the pyrolysis rate of each experiment shows sustained 
fluctuations. 


The next two figures, Figures 5.2 and 5.3, show the upper and lower layer temperatures for the 
fire room. The figures display thermocouples CT4543 and CT4514 which are respectively 
located at the ceiling and floor of the fire room approximately two meters from the wood crib 
platform. 
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Figure 5.2: Fire Room Upper Layer Temperature (CT4543) 
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Figure 5.3: Fire Room Lower Layer Temperature (CT4514) 
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The following conclusions can be drawn for the upper layer temperature characteristics from 
Bigure:s2; 


e All three tests depict an identical temperature behavior during the first six minutes after 
ignition. This most likely results from the fire at this stage being consisting solely of wood 
crib number 2 and the alcohol source used for ignition. | 

e While experiment T51.16 reaches a plateau at 300 °C, lasting for about 2 minutes, 
temperatures for T51.17 and T51.18 continuously increase with T51.18 showing slightly 
higher temperatures. 

¢ At 10 min., all three tests undergo a very steep temperature gradient (flashover), with T51.16 
depicting a slight time delay of 1 min. because of the preceding temperature plateau. 

e At 15 min., test T51.18 has a large drop in the temperature gradient corresponding to the 
underventilation period shown in Figure 5.9. This new gradient lasts until the peak fire 
temperature is reached. 

¢ Tests T51.16 and T51.17 both reach peak temperatures at 20 min. Test T51.18 reaches its 
peak temperature at 28 min. 

¢ From 30 min. till the end of the tests, all the experiments show the same cooldown behavior. 


From Figure 5.3, depicting the lower layer fire room temperatures for the three tests, the 
following observations hold: 


¢ All three tests show identical temperature behavior, a benign and slow increase over the first 
6 min. 
¢ The same characteristics can be identified for the lower layer temperatures as for the upper 
layer temperatures beyond 6 min. 
¢ Tests T51.17 and T51.18 show identical behavior; their temperature gradients in the lower 
layer are higher than in the upper layer, after flashover occurs, while the temperature 
gradient for T51.16 seemingly less than in the upper layer and obviously less than the other 
two tests. 
¢ The maximum temperatures are reached at somewhat later time than for the upper layer. 
¢ The following peak temperatures are reached: 
T5116:280"°C at22 min: 
T51.17: 480 °C at 24 min. 
TS118: 950" Crar29 mint 
e The temperature trace of CT4514 for test T51.16 and to some extent for test T51.17 indicates 
distinct fluctuations over the time period of 16 to 35 min which is indicative of a mixing 
layer. 


The next two figures on the following page, Figures 5.4 and 5.5 depict temperatures in the 
hallway connecting the fire room to the curtained area. The sensors are located at the midpoint 
of the hallway both in terms of length and width. They are at elevations approximately 0.5 m 
below the ceiling and 0.5 m above the floor. 
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Figure 5.4: Hallway Upper Layer Temp. (CT4673) 
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Figure 5.5: Hallway Lower Layer Temp. (CT4663) 
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The following items are observed in these figures: 


e As in the fire room, the hallway temperatures in both the upper layer and the lower layer 
show a continuing increase with the number of wood cribs, but the maximum temperatures 
are even more distinctly separated. 

¢ While T51.18 reaches a maximum of 680 °C at 23 minutes, experiment T51.17 reaches only 
430 °C and T51.16 reaches 310 °C at the same time, e.g. times of peak temperatures of the 
upper gas layer in the hallway at that position coincide for all three tests. 

¢ The percent change in temperatures between tests T51.16 and T51.17 was 60% for the lower 
layer and 40% for the upper layer. For tests T51.17 and T51.18 these values were 80% and 
60% respectively. The change in percentage increase between T51.16/T51.17 and 
T51.17/T51/18 is the same for both layers, 20%. 

¢ The hallway temperatures have nearly identical trends as a function of time with the 
exception of overall magnitude. 

e The significantly cooler temperatures in the lower layer indicate the presence of a large flow 
of fresh air under the maintenance hatch curtain to the fire room. 


The next three figures, Figures 5.6 through 5.8, show temperatures in the main staircase 
maintenance hatch on the 1.400 level and the 1.700 level and in the dome. Thermocouple 
CT4671 is located approximately 0.5 m below ceiling of the 1.400 level, CT7703 is located at 
the middle elevation of the 1.700 level, and CT0423 is located 3 m above the dome operating 
deck. The comparisons between these three figures as well as Figures 5.2 and 5.4 lead to the 
following observations: 


¢ The upper gas layer temperatures for all three experiments drastically decrease along the 
flowpath from the fire room to the level 1.400 main staircase maintenance hatch (CT4671), 
such that for experiment T51.18, 11 cribs, the hot layer temperature drops from 1000 °C to 
350 °C and similarly for the other two tests. 

¢ The maximum temperatures are reached at the same time, 24 min, for all three experiments 
for the upper layer at the 1.400 level maintenance hatch. 

¢ While a pronounced difference of 140 °C exists between the peak temperature between 
T51.17 and T51.18 much less of a difference, 50 °C, remains between the former and 
T51.16. In fact, both of these experiments indicate about the same cooldown behavior. 

¢ The hot plume rising up the maintenance hatch cools quickly between the 1.400 level and the 
1.700 level which indicates a large entrainment of air from the upper levels into the plume, 
by both entrainment via counter-current flow of cold air descending through the hatches and 
by entrainment from the surrounding compartments. As a result, the differences between the 
peak temperatures decrease and the overall appearances of the temperature traces become 
much more similar, Figure 5.7. Heat transfer to the structures also contributes toward 
cooldown. | 

¢ The temperature data for the plume shows oscillations of up to 20 °C at the entrance to the 
1.500 level, Figure 5.6, demonstrating the turbulent nature of the plume. 

¢ Between the 1.700 level and the dome the entrainment reduces due to limitation by the flow 
area restrictions in the maintenance hatch. Consequently, the cooldown of the plume is only 
minimal, a few °C, but the ranges of temperature increases are further broadened over time 
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while the overall shapes of the temperature curves become more similar with even further 
reduction of the original differences between T51.17 and T51.18, Figure 5.8. 
Although greatly dampened by entrainment and heat transfer, the plumes entering the dome 


are distinctly identifiable for all three tests despite the long complex flow path originating in 
the fire room , 30 m below. 


400 T5iniG 
5 cribs 
350 
Slee 
mun 300 7 cribs 
O 
© 250 11 cribs 
© 200 
® 
(on 
£6150 
o 


100 
50 


rom O10 20° 307-40 50° 60 70 80790 
Time (min) 


Figure 5.6: Level 1.400 Main Staircase Maintenance Hatch Temp. (CT4671) 
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Figure 5.7: Level 1.700 Main Staircase Maintenance Hatch Temp. (CT7703) 
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Figure 5.8: Dome Temp. (CT0423) 
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Figures 5.9 through 5.11 present the oxygen, carbon monoxide, and carbon dioxide 
concentrations in the upper layer of the fire room doorway. The following observations can be 
drawn: 


¢ A steady drop in the minimum oxygen concentration is seen with an increase in the number 
of wood cribs. 

¢ Test T51.18 reaches 0% oxygen concentration for a period of 15 minutes which indicates 
this fire is clearly underventilated. 

¢ Tests T51.16 and T51.17 generate only minute amounts of CO while test T51.18 produces a 
significant amount of CO. This also indicates that T51.18 quickly became underventilated 
while the other two wood crib tests were able to maintain a sufficient airflow into the fire 
room from the hallway to completely combust the wood. 

¢ CO, concentrations echo the reverse of the trends seen in the O2 concentrations. A steady 
increase with the number of wood cribs and a plateau in the values for test T51.18 that 
corresponds to the period of complete oxygen depletion in the fire room. 

¢ Many of the characteristic details seen in the temperature histories are reflected in the gas 
concentrations. Thus, they are confirmed independently by means of diverse measurement 
principles. 

¢ COconcentrations in tests T51.16 and T51.17 peaked at 70 ppm and 100 ppm respectively, 
well below the minimum 30 minute exposure concentration for lethality which is 4000 ppm 
[35]. Test T51.18, however, reached a CO concentration of 30000 ppm or well above the 
minimum lethal concentrations. 
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Figure 5.9: Fire Room Doorway Upper Layer O, Concentration (CG4642) 
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Figure 5.10 Fire Room Doorway Upper Layer CO Concentration (CG4641) 
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Figure 5.11: Fire Room Doorway Upper Layer CO, Concentration (CG4643) 
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The following two figures, Figures 5.12 through 5.14, show the measured velocities in the upper 
layer of the fire room doorway, the lower layer of the hallway near the curtained area, and in the 
main staircase maintenance hatch at the 1.700 level. 


¢ All of the experiments reach a maximum velocity of about 6 m/s in the upper doorway 
region and maintaining that level over a time span of between 14 min. , T51.16, and 27 min., 
Polts: 

¢ Surprisingly, the highest measured velocity is for test T51.16, the lowest fire power of the 
wood cribs. This is probably due to the location of the peak of the velocity profile in the 
doorway. The velocity probe is located 0.17 m below the top of the doorway. For test 
T51.16, which had a smaller hot layer, the peak of the velocity profile was probably closer to 
the probe than during the other two tests. 

¢ Measured velocities exiting the fire room are very similar for each of the tests. 

¢ The change in the upper layer temperature gradient in the fire room during test T51.18 is 
also seen in the upper layer velocity for that test. 

¢ Since the hallway sensor for test T51.16 failed, only tests T51.17 and T51.18 can be 
compared. The peak velocity in the hallway increased from 3.5 m/s to 5.1 m/s for these two 
tests, an increase of 44%. This compares well to the increase in the number of wood cribs 
from 7 cribs to 11 cribs which is an increase of 57%. 

¢ A steady increase in the velocities measured in the maintenance hatch is seen for each 
increase in the number of wood cribs. 

e The fire room and hallway velocities show some high frequency fluctuations which are 
indicative of turbulent flow. In the maintenance hatch, these high frequency fluctuations are 
superimposed on top of low frequency fluctuations. These low frequency fluctuations which 
can take minutes to occur suggest that the plume rotates its position about in the hatch. This 
behavior is considerably more complex than the simple static position plume models used in 
fire codes. 

¢ All three experiments start with an initially existing velocity field in the containment with 
fire room and hallway velocities of about 0.5 m/s. T51.16 at the 1.700 level maintenance 
hatch depicts stagnation. 

¢ The hallway lower velocities of 3.5 m/s and 5.1 m/s for experiments T51.17 and T51.18 
indicate a strong counter-current flow of fresh air with nearly as high of velocities as seen in 
the upper layer. 

¢ At the 1.700 level, peak velocities range between 3 m/s for T51.16 and 3.9 m/s for T51.18, 
e.g. a substantially high value. 

¢ At 90 min., well after the end of combustion, velocities in the flow path range well above | 
m/s and 2.5 m/s depending on the proximity to the fire room, where they are higher. 

¢ Also, details in the velocity traces independently confirm major characteristics observed in 
the temperature histories in the fire room. 
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Figure 5.12: Fire Room Doorway Upper Layer Velocity (CV4640) 
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Figure 5.13: Hallway Lower Layer Velocity (CF4612) 
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Figure 5.14: Level 1.700 Main Stairchase Maintenance Hatch Velocity (CF7701) 
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Figures 5.15 and 5.16 on the following page present the measured CO, concentrations in the 
main staircase maintenance hatch at the 1.600 level and in the dome. 


The CO, trends for both locations are similar in both magnitude and shape. The CO 
concentration in the dome shows a spreading of the peak at 25 minutes compared to the 
hatch, but this is to be expected given the large free air volume present. 

The final, steady state concentrations at both locations are identical indicating that the flows 
induced in the containment lead to a total mixing of the containment atmosphere, 
homogenizing the CO, concentrations. 

The percent increase in the steady state CO, concentration from test T51.16 to test T51.17 
was 35%. This matches closely the increase in the amount of fuel, 5 cribs to 7 cribs or 40%. 
Given that both of these tests were well ventilated, this serves to qualitatively confirm the 
sensor readings. The percent increase between test T51.17 and test T51.18 is slightly lower 
than the increase in the number of wood cribs, a 42% CQ» increase compared to a 57% fuel 
increase; however, this also acts to confirm the sensor readings and the general 
underventilated nature of test T51.18. 


The final two figures of this subsection, Figures 5.17 and 5.18, show the extinction coefficients 
in the main staircase maintenance hatch on the 1.700 level and in the dome. 


There is a fairly small increase in the extinction coefficients between test T51.16 and T51.17, 
since both of these tests were rather well ventilated. Contrary, test T51.18 shows a very 
large increase that clearly demonstrates the tremendous amounts of smoke produced during 
this test. 

Smoke density drops by nearly a factor of two from its maximum value as the containment 
atmosphere continues to mix, cooldown, and particulate matter settles out onto the building 
surfaces. 

Assuming a visibility limit of 0.5 for trained rescue personnel familiar with the facility [8], 
test T51.18 very quickly reached a point where rescue operations within the facility would 
have become impossible. 

For untrained personnel, assuming a visibility limit of 0.15 [8], it would have been 
impossible to successfully evacuate for the entire T51.17 and the T51.18 tests and during the 
peak combustion portion of the T51.16 test. 
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Figure 5.15: Level 1.600 Main Staircase Maintenance Hatch CO, Conc. (CG1166) 
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Figure 5.16: Dome CO, Concentration (CG1104) 
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Figure 5.17: Level 1.600 Main Staircase Maint. Hatch Extinction Coeff. (CG2166) 
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Figure 5.18: Dome Extinction Coefficient (CG2104) 
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5.2 Representative Experimental Results 


This subsection gives the value of various parameters for each of the wood crib tests at 16 
minutes after the start of the fire. This time represents the peak pyrolysis rate for all three tests. 
The figures are plotted as parameter vs. number of wood cribs combusted. This allows for an 
easy comparison of the parameters on the basis of available fuel. Figures 5.19 and 5.20 show 
the upper and lower layer temperatures in the fire room. The following observations are made 
from these figures: 


¢ Upper layer temperatures increase between tests T51.16 and T51.17. Little change in layer 
temperatures is seen between tests T51.17 and T51.18. 

¢ Lower layer temperatures increase between test T51.16 and T51.17. During test T51.18 the 
lower layer temperatures do not increase further. 

* Both layers have a non-uniform temperature distribution clearly showing the effects of the 
doorway, the location of the fire, and a cold region along the wall furthest from the wood 
cribs. 

¢ Actually, CT4512 reaches the highest temperatures recorded for T51.16 with five wood cribs 
in the lower layer. 

¢ The maximum temperature difference between the upper and lower layers for all three tests 
is identical, 700 °C. Since this difference remains constant, rather then increasing, with the 
number of would cribs it indicates that the tests become increasingly ventilation controlled as 
the number of wood cribs increases. 

¢ Given the large overlap of temperatures between the two layers, and the temperature 
distribution itself, the concept of layers seems to be questionable at this point in the tests. 


Figures 5.21 and 5.22 compare temperatures along the upper and lower layers of the hallway 
leading from the fire room to the maintenance hatch. Most of the observations made above for 
the fire room also apply to the hallway. In addition the following conclusions can be drawn: 


¢ Both layers show nearly linear increases in temperatures with the number of wood cribs. 

¢ The upper layer plume cools significantly as it moves down the hallway; however, during 
test T51.18 the cooling of the plume does not begin until a point halfway down the hallway 
as seen by the nearly identical temperatures for CT4674 and CT4673. Due to the 
underventilated nature of T51.18 burning may have occurred in the hallway beyond the fire 
room doorway. This would prevent the cooling of the plume in the hallway. 

¢ The air entering the hallway’s lower layer is just above the ambient temperature. This cold 
air is coming underneath the curtain surrounding the maintenance hatch. Its near ambient 
temperature indicates that the hot plume from the fire room does not mix with the cold 
returning air in the curtained area. 

¢ The air exiting the hallway’s lower layer into the fire room is significantly above the ambient 
temperature indicating a large degree of thermal mixing between the hot and cold gas 
streams occurs in the hallway. 
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Figure 5.19: Fire Room Upper Layer Temperatures at 16 Minutes 
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Figure 5.20: Fire Room Lower Layer Temperatures at 16 Minutes 
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Figure 5.21: Hallway Upper Layer Temperatures at 16 Minutes 
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Figure 5.22: Hallway Lower Layer Temperatures at 16 Minutes 
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The next two figures, Figures 5.23 and 5.24, show temperatures in the main staircase 
maintenance hatch and in the dome. The following observations are made: 


¢ CT5303 which lies in the open area of the maintenance hatch on the 1.500 level is seemingly 
located on the edge of the plume because its temperatures are lower than the thermocouple 
above it, CT6609 for all three tests. 

¢ The plume entrains significantly and cools between the 1.500 level and the 1.600 level; 
however, it cools only slightly between the 1.600 level and the 1.700 level. At some point 
on the 1.600 level the plume becomes constrained due to the flow restrictions created by the 
walls and the limited opening of the maintenance hatches, it can partly spread over to the 
staircase. 

¢ Within the dome the plume cools to ambient temperatures as it reaches the spiral staircase. 
The large free air volume in the dome allows for significant entrainment into the remaining 
buoyant plume. 

¢ All the temperatures show an increase with the number of wood cribs; however, this increase 
is not linear due to the underventilated nature of test T51.18. 


Figures 5.25 and 5.26 display the O» and CO, concentrations measured in the upper layer of the 
fire room doorway. As expected, the data indicate that as fire power increases the hot gasses 
leaving the fire room become increasingly depleted in O» and increasingly enriched in CO) at the 
same time. The underventilated nature of test T51.18 is clearly indicated by both the 0% O, 
concentration and by the nonlinear increase in CO» concentrations as the number of wood cribs 
increases. 


The final figures in this section, Figures 5.27 and 5.28, show CO, concentrations at various axial 
levels in main and spiral staircase maintenance hatches. The following is noted: 


¢ CG1146 in the main staircase is located in the lower layer of the curtained area. It shows 
nonzero CO; levels in the airflow beneath the curtain. This results from the global 
circulation of combustion products throughout the HDR facility as a result of the vertical 
shafts formed by the staircases and maintenance hatches. At 16 min. combustion products 
have been transported up to the dome via the main staircase flowpaths and back down to the 
lower levels by the spiral staircase flowpaths. 

* CG1153 does not appear to be within the plume as it shows much lower CO, levels than 
other instruments in the main staircase region. 

¢ As already seen in the main staircase temperatures, Figure 5.22, also confirms that the CO, 
levels decrease insignificantly in the upper levels of the facility just below the dome. This 
further supports the conclusion that the plume has become limited by flow restrictions and 
bounding walls and that the entrainment is small little air in that region. 

¢ A continuous drop in the CO, concentration is achieved as the plume circulates through the 
HDR facility. 

¢ Dilution of the CO, concentration between 10.2 m and 31.0 m along the plume is 0.14 vol % 
to 0.17 vol %. 

¢ CO, concentration decreases by a factor of two to three between the main and the spiral 
staircases due to mixing and entrainment in the dome. 
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CO), concentration increases in the spiral staircase in the direction of lower axial positions, 
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Figure 5.23: Main Staircase Hatch Temperatures at 16 Minutes 
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Figure 5.24: Dome Temperatures at 16 Minutes 
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Figure 5.25: Doorway Upper Layer O, Concentration at 16 Minutes 
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Figure 5.26: Doorway Upper Layer CO, Concentration at 16 Minutes 
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Figure 5.27: Main Staircase Maint . Hatch CO, Concentrations at 16 Minutes 
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Figure 5.28: Spiral Staircase Maint. Hatch CO, Concentrations at 16 Minutes 
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6 POTENTIAL OF EXPERIMENTAL DATA FOR CODE VALIDATION 


One of the primary purposes of the HDR fire experiments was to create a database of 
experimental data for use in code validation and model development. This section will discuss 
aspects of the T51 wood crib fire tests and the related data that can be used for code validation. 
The tests will be discussed in terms of zone models, containment system codes and field models. 


6.1 Zone Models 


Zone model fire codes, such as CFAST [29,30] and MRFC, operate by assuming that in a fire 
situation every room in a building can be represented by two layers: a hot layer containing the 
combustion products from the fire and a cold layer which is oxygen rich. A number of elements 
of the instrumentation plan for the T51 wood crib tests were established for the purpose of 
collecting data for the evaluation of zone model codes. Some of these elements are discussed in 
this subsection. 


6.1.1 Layer Height 


The key parameter calculated by a zone model code is the layer height. The instrumentation 
plan for the wood crib tests was such that layer heights are only obtainable in the fire room and 
the fire room doorway. The fire room was instrumented with an array of thermocouples that had 
four axial levels and the doorway with a rake of thermocouples having five axial levels. The 
remainder of the facility lacked vertical rakes of thermocouples. 


6.1.2 Layer Temperatures 


For the three wood crib fires sufficient thermocouples are available to evaluate the layer 
temperatures in the fire room, the fire room doorway, and the hallway. The fire room 
thermocouples were discussed above. In the hallway instrumentation consisted of three sets of 
thermocouples along its length. Two of these were triplets of thermocouples with one near the 
floor of the hallway, one near the hallway’s middle elevation, and one near the hallway’s ceiling. 
The third set was a pair of thermocouples with one at the floor and the other in the middle of the 
hallway. Thus, the temperature for each layer in the hallway can be obtained an average of the 
temperatures for the entire length of the hallway. 


6.1.3. Mass Flow Rate 


In those locations were layer height and velocity information exist or where a velocity sensor 1s 
located in a maintenance hatch, mass flow rates between compartments or between levels can be 
determined. In the case of horizontal flow the doorway dimensions and layer height information 
is used to determine the flow area of the layer. This along with the ideal gas law and layer 
temperature then yields the mass flow rate. For vertical flow, assuming that the plume occupies 
the whole hatch can also yield a rough estimate of the mass flow estimate. 
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During the wood crib fire tests T51.1 experiments only two velocity sensors yield horizontal 
mass flow data. These were located in the upper layer of the fire room doorway and in the lower 
layer of the hallway. 


A number of sensors were positioned in the hatches for the wood crib fire tests. Mass flow data 
is therefore available for the main staircase on three levels: the 1.500 level, the 1.700 level, and 
the 1.800 level. Data are also available for flow down the spiral staircase hatch on the 1.700 
level. 


6.1.4 Optical Smoke Densities (Extinction Coefficients) 


The wood crib fire tests produced varying quantities of smoke which were transported 
throughout the HDR facility. Test T51.18 produced large quantities of smoke which heavily 
coated the surfaces in the facility. A number of sensors were placed in the facility to measure 
the optical smoke density. The sensors were in three locations: a vertical rake in the hall 
connecting the spiral and main staircases on the 1.600 level, in the main staircase hatch on the 
1.600 level, and in the main staircase hatch in the dome. 


6.2 Containment System Codes 


Containment system codes, such as GOTHIC [31] and RALOC, owe there origin to the nuclear 
power industry. The need to evaluate the effects of loss-of-coolant accident scenarios in a 
containment building requires thermal hydraulic codes that accurately model the two-phase, 
thermal-hydraulic response of a large building to a source of energy, mass, and momentum, e.g. 
a break in a reactor coolant pipe. Containment system codes commonly use the lumped 
parameter method, that is a compartment is considered to be a single point whose properties 
represent the volume-averaged properties of the compartments. By modifying the source to be 
combustion gases and radiant heat rather than steam and water, containment system codes can be 
applied to computing the effects of fires on large structures. These codes usually model all vents 
between compartments, heat transfer to structures, sprays, ventilation systems, etc. In the case 
of GOTHIC, the discretization options also include a combination of lumped and distributed 
parameter nodalizations (1-D, 2-D, and 3-D). 


For the purpose of fire modeling some adjustments to the typical lumped parameter approach 
must be made. In order to appropriately generate the buoyant driving force for fire drive flow, 
the fire compartment and any immediately adjoining rooms cannot be modeled as lumped 
volumes. Rather each room must be subdivided into a network of lumped volumes to allow for 
thermal stratification in these compartments. 
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6.2.1 Compartment Temperatures 


As many regions outside the fire compartments and adjoining compartments are modeled as 
lumped volumes, this results in more useable temperature information in the T51 data set than 
for the zone models. For all of the T51 tests, temperatures can be obtained for the fire room, the 
doorway, the hallway, the curtained area, the 1.400 level outside the curtain, at each level for 
both the ascending and descending flow through the maintenance hatches, the dome, and the 
connecting hallway between the spiral and main staircases. These regions represent the entire 
circulation loop induced by the wood crib fires. 


6.2.2 Compartment Mass Flows 


For the most part, system codes rely on single or multiple openings to connect compartments. 
Therefore, the wood crib tests are not very useful for system codes for comparing mass flows 
due to the lack of velocity sensors. The vertical flows up the hatches are suitable for comparison 
with system codes. , 


6.2.3 Wall Heat Conduction 


The wood crib tests featured heat transfer blocks. The data from these sensors have been used to 
determine the time-dependent heat transfer into the HDR surfaces. As the storage and release of 
energy into structures is very important in nuclear accident analysis, system codes tend to 
contain robust algorithms for heat conduction into layered structures, e.g the coated and painted 
wall surfaces. The measurement blocks used in the HDR were designed for the purpose of 
evaluating these algorithms as well as obtaining the steady-state and the transient behavior of the 
heat transfer coefficient to compare with common correlations. 


6.2.4 Combustion Models 


The nuclear containment codes were originally designed for the purpose of calculating a 
containment building’s response to a loss of coolant accident in a nuclear power plant. 
Therefore, the energy, mass, and momentum sources that were originally coded were 
steam-water mixtures and hydrogen from metal-water reactions. For these codes to be used to 
model the effects of fires in containment buildings, additional models relating to the calculation 
of combustion physics were added. Unlike the models relating to energy and mass inputs from 
pipe breaks, these new combustion models have not been thoroughly tested. Therefore the 
pyrolysis and combustion product measurements made during these tests are invaluable to 
evaluating the effectiveness of these models. 
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6.3 Field Models 


Field models, such as FLUENT [36], FLOW-3D, and NIST-LES [37], operate by solving a 
discretized form of the three-dimensional equations for mass, momentum, and energy 
conservation. For most real structures, accurate resolution of the velocity, temperature, and 
species field for a fire require a large number of computational nodes. Therefore, use of field 
models is typically restricted to smaller subsets of a larger structure to reduce the computational 
resource requirements. 


The wood crib fire tests were sparsely instrumented for the needs of a field model; however, the 
fire room instrumentation is sufficient for some field model testing. Since the wood cribs were 
located in front of the fire room doorway, the region of the fire room without the wood crib was 
a cooler region of the fire room. A field model code should be able to reproduce that local 
phenomena. Also, the field model should be able to reproduce the temperature distribution seen 
as a result of the fire itself and due to the inrush of cooler air from the hallway. This information 
along with velocity and gas sensors in the doorway can be used to qualitatively and for some 
aspects quantitatively verify a field model of the fire room. The adequacy of this approach may 
depend on the accuracy of the boundary conditions to be applied to the model based on rather 
coarse instrumentation. 
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7 T51 WOOD CRIB CFAST MODEL 


The T51 wood crib fire tests were modeled using a modified version of the C model used for the 
TS51 gas fire tests [34]. This section will discuss the runtime and environmental parameters, the 
combustion properties for the fire, the room surface properties, and the geometric model 
developed for the CFAST computations [29,30]. A full description of the development path for 
this CFAST model can be found in Volume 2 of the report [34]. It is noteworthy to mention that 
the most significant change to the gas fire model, other than the fuel, is that the upgraded version 
of CFAST used for this report, v3.1.1, was capable of modeling the HDR as an enclosed facility 
unlike the gas fire model which unrealistically required a vent open to the outside for the 
duration of the run. 


eh. Environment and Runtime Control 


The T51 wood crib tests had a fire duration of approximately 45 minutes. Data collection during 
the tests began 10 minutes prior to the fire’s start and continued for a period of 40 to 60 minutes 
after the end of the fire. No preconditioning of the containment building was performed before 
the start of any tests, 1.e. the containment was at ambient conditions at the start of each test in the 
series. The CFAST model reproduced these conditions with a 100 s pretest period with no fire 
to generate a baseline followed by the appropriate fire conditions for a duration of 90 minutes 
from the ignition of the wood cribs. Table 7.1 below shows the CFAST input cards for 
environment and runtime control. 


Table 7.1: Environment and Runtime Control Input Cards for CFAST 


Simulation Time (s) 100 s pre-fire + 90 min post-fire 
Print Interval (s) 

History Interval (s) 

Display Interval (s) ae graphics 


Cop y Count_ sraphics to cop 


Ambient Pressure (Pa) l ho 
Station Elevation (m) ____| Set hallway floor as 0 elevation 


1p} Combustion Model 


Each of the tests in the T51 wood crib series used a different number of identically constructed 
wood cribs set on a platform in the fire room, see Figure 2.2 on page 2-5. This platform was 
placed on top of three digital scales from which the pyrolysis rate was determined. These 
measurements form the basis for the combustion parameters used for the CFAST models. 


Wood has an approximate empirical chemical formula of C 0.9; H240. The heat of combustion of 
the TYPE OF WOOD used in the T51 wood crib tests is estimated to be 12600 kJ/kg. This 
estimate is based on the maximum fire powers developed during the wood crib tests as reported 
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in [10]. However, this empirical formula cannot be directly used in CFAST v3.1.1. Specifying 
oxygen as part of a fuel description in CFAST results in the formation of unrealistic combustion 
products. This is discussed in Appendix B. Therefore, the oxygen in the wood was ignored and 
the pyrolysis rate and heat of combustion adjusted as follows to conserve species mass and 
energy input for the models: 


¢ Wood is 38.3% C, 21.1% H, and 53.6% O. Therefore, when the oxygen is considered absent 
53.6% of the experimentally pyrolyzed mass must be removed in order to keep constant the 
mass of carbon and hydrogen available for combustion. The experimentally measured 
pyrolysis rates were reduced by 53.6% before being input into CFAST. 

¢ Ifthe mass of fuel is reduced by 53.6%, the heat release rate will also be reduced if the heat 
of combustion is not adjusted. To preserve total energy release, the heat of combustion was 
increased by 116% to 27200 kJ/kg. 


Testing on these changes was performed using a single room model to verify that the appropriate 
total heat release and CO, was achieved. 


The final parameters required for specifying the fuel was to provide the CO and soot production 
values. In reality these parameters are quite complex for a wood fire, and given the available 
data for the HDR, they were not known. Therefore, the CFAST preprocessor was used to define 
a single compartment fire using plywood as a fuel. The preprocessor determined values were 
assumed to be reasonably representative of wood crib fires and thus subsequently used for input 
into CFAST. 


With the combustion product information and the basic fuel information, the CFAST input cards 
relating.to combustion were completed. The applicable cards and their values are shown in 
Table 7.2 for test T51.16 which was a test with 5 wood cribs. For the other tests modeled, 
modified cards include FAREA (changed to reflect the number of wood cribs used), FMASS (to 
reflect the measured pyrolysis rate for each test), and FQDOT (to reflect the current test’s fire 
power). These fire power related changes are documented in the input decks presented in 
Appendix A. The remaining combustion parameters were not changed between the tests 
modeled. 
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7.3 Surface Properties 


To model the HDR with CFAST requires generating material properties for the HDR’s 
construction materials. For the T51 wood crib tests there are three materials which must be 
defined which are the HDR structural concrete, Ytong fire brick, and Alsiflex fireproof, glass 
fiber matting. A CFAST material library was created for these materials and the surfaces which 
were constructed from them. Table 7.3 gives the thermophysical material properties which were 
obtained from [20]. 


Table 7.3: Thermophysical Properties for HDR Construction Materials 
Density 


Banas K) are rs = 
i 0p 
mea aa re 
ae Pe 
| Alsiflex Mats _| 


' Taken as firebrick from wea eae material ioraryee tore library 
* Taken from Alsiflex properties in [35] 


The above thermophysical properties were used to create a seven material library for use with 
CFAST. The seven materials included five, single layer materials and two, two layer materials. 
The materials in the order listed below are 3 cm thick Alsiflex matting, 100 cm of concrete, 50 
cm of concrete, 10 cm of Ytong firebrick, 25 cm of Ytong firebrick, the Alsiflex matting and 
Ytong firebrick ceiling for the fire room, and the Ytong firebrick and concrete floor of the fire 
room. The material library is shown below. For single layer materials the format is material 
name, conductivity, density, heat capacity, thickness, emissivity, and seven parameters for HCl 
production. For multiple layer materials the format is the same with a ‘/ denoting values for 
each layer. 


ALSTELEX 50 al ToL000P LS 0R0 10 3202970 "OE OF US0-80 D 

CONCRI00 S42 1 “S79 2225. 102 820 0480 040. O20 

CONCRO5O jv2. 1 87S, 2225 Wiss OF on 0 SO 30 10? 02070 

YTONGLOO “O'24" 950" 340, Oa 0S OO (OTe Ceti e meni 

YTONGZ50. 024. 950) (340° 0225) Use m0 0G Ome eG 

FIRECHIL, 0.17 /D0e24 2000/950 1307340400370 25505970 0 OOm0 0: 0 
IRE FUR 0.24/2.1 950/679 3400/2225 0.25/1 Osep Ore U2020806e8 


7.4 Compartments and Compartment Interconnections 


The final portion of the CFAST model is the compartment definitions and the compartment 

interconnections. As the HDR facility contains 9 levels and over 60 compartments with complex 
interconnections, it is not feasible to model it explicitly with CFAST. Therefore, simplifications 
were made to model the HDR facility with CFAST. These simplifications reduced the HDR to a 
three level structure composed of the fire compartments on the 1.400 level, the 1.600 level, and a 
lumped level consisting of all building elevations above the 1.600 level. The development of the 
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CFAST model benefited greatly from the experience gained simulating the gas fire tests [34] 
plus some progress made in version 3.1.1 relating to completely enclosed fires. 


Figure 7.1 shows the 1.400 level. From this figure it is obvious that the important compartments 
were the fire room and its doorway, the hallway, and the curtained area as these compartments 
defined the location of the fire and the flowpath taken by the combustion products. The 
curtained area had a 0.5 m gap at the floor for the return of cold air. This air could come from 
either room 1.402 or the main staircase, which in itself was connected to room 1.402. Room 
1.402 had a fairly narrow connection to the rest of the containment building, so this restriction 
should be accounted for to create a flowpath resistance for returning cold air. Therefore, the 
1.400 level of the model consisted of five compartments which were the fire room and its 
doorway, the hallway, the curtained area, the main staircase, and room 1.402. 


From the dimensions given in Figure 7.2, the volume of the fire room and doorway was 
calculated to be 29.6 m* with the doorway having 3 m’ of that volume. Since the doorway was 
such a small volume in comparison to the fire room and the hallway, it was assumed to be part 
of the fire room. Therefore, as the fire room was 2.75 m in height, the floor area of the fire 
room is 10.77 m?. If the linear distance from the fire room’s back wall to the doorway opening 
into the hallway was preserved, as this is represents the room dimension about the wood crib 
platform, the floor dimensions of the fire room became 3.65 m x 2.95 m. The surfaces of the 
fire room consisted of a ceiling lined with 0.25 m of Ytong firebrick covered by 0.03 m of 
Alsiflex mats, the walls of the fire room lined with 0.25 m of Ytong firebrick, and the fire room 
floor which was 0.25 m of Ytong firebrick on top of a 1 m thick concrete floor. The fire room 
was connected to the hallway by a 1.01 m x 1.975 m doorway. 


0.250 m 


= 
8 : [6] F 
Gas Burners (6) ee i e [9] 5 
8 So 4 
180° 3S Ee 9 
Figure 7.1: Level 1.400 Figure 7.2: Fire Room and Doorway 
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Figure 7.3 shows the hallway that connected the fire room to the maintenance hatch. Reference 
[9] gives the volume of the hallway as 22.15 m° with a ceiling height of 2.485 m. The hallway 
had a variable width and a curved length which did not readily translate into a rectangular 
parallelepiped for modeling purposes with CFAST. The dimensions of significance for the 
hallway were the height and the width of the hallway where it connected to the maintenance 
hatch; this was also the widest cross section of the hallway. This resulted in the hallway having 
floor dimensions of 1.8 m x 4.95 m. The hallway surfaces were 0.10 m of Ytong firebrick for 
the walls and ceiling and 0.50 m of concrete for the floor. The hallway was connected to the 
curtained area by a 1.8 m x 2.75 m opening. 


Figure 7.4 shows the curtained area beneath the maintenance hatch. The hatch itself was a 
rectangular opening 4.54 m? in area. Using the hatch dimensions shown in the figure and 
assuming the figure was to scale, the curtained area floor dimensions were found to be 4.30 m x 
2.75 m. The height of this region was 4.60 m. All surfaces were assumed to be 0.50 m thick 
concrete. The hatch was connected to the outside by a 4.54 m’ rectangular opening. From 
Figure 7.4, the flow path underneath the curtain to room 1.402 and the main staircase was 
estimated to be 2.0 m x 0.5 m and 4.3 m x 0.5 m, respectively. 


lo) 


= 


@ 
Doorway 


Hallway 


Figure 7.3: Hallway Figure 7.4: Maintenance Hatch and Curtain 


The main staircase had the same floor and ceiling elevations as the hatch area, 4.60 m in height. 
The floor area of the main staircase was estimated to be 11.78 m’, 4.33 m x 2.72 m, by using 
Figure 7.1 and the HDR steel shell diameter of 20 m. All surfaces were again assumed to be 
0.50 m thick concrete. The vertical flow path up the main staircase was 5.75 m*. The 
connection to room 1.402 was estimated from Figure 7.1 to be 2.7 m x 3.1 m. 


Since experimental data from the T51 test series indicated that the 1.500 level was not involved 
in the global flows during the tests it was not modeled in CFAST. Therefore, the 1.400 level 
vertical connections from the hatch and staircase had their ceiling heights extended to connect 
them directly to the 1.600 level. This resulted in the compartment height for these two 
compartments to be implemented as 11.10 m rather than 4.60 m. 
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Room 1.402 also had its floor area dimensions estimated from Figure 7.1. The floor dimensions 
were estimated to be 1.80 m x 6.50 m with a height of 3.50 m. All surfaces were 0.50 m thick 
concrete. The connection from room 1.402 to the outside is given in [1] as 0.8 mx 1.5 m. 


On the 1.600 level, there were three significant regions, as shown in Figure 7.5. The vertical 
flow paths created by the 1.400 level main staircase and maintenance hatch must be continued. 
The portion of the 1.600 level connecting the two staircase regions needed to be included. 
Finally, the spiral staircase and its maintenance hatch also needed to be included. Since room 
1.603 is an isolated compartment that is connected solely to the spiral staircase it can be 
considered part of the staircase region. Therefore, the 1.600 level portion of the model 
contained four compartments: the curtained area on the 1.600 level, the main staircase on the 
1.600 level , the spiral staircase and room 1.603, and the remainder of the 1.600 level. 


The main staircase and hatch compartments were given the same floor dimensions on the 1.600 
level as they were on the 1.400 level. The height of these compartments was set to 5.05 m 
which places the ceiling of these compartments at the elevation of the 1.700 level’s floor. All 
surfaces were set to 0.50 m thick concrete. The maintenance hatch area leading to the outside 
remained unchanged from the A model, but the main staircase vent was set to a 6.97 m° 
rectangle. 


The volume of the room 1.603 and spiral staircase compartment was set to equal the sum of the 
compartment volumes for 1.603 and 1.611 which is 472 m’*. A height of 5.05 m was also used 
for this compartment. The dimensions of the floor were set to a square as there was no dominant 
linear dimension of interest for this compartment. The floor area was therefore set to 9.67 m x 
9.67 m and all room surfaces were 0.50 m thick concrete. The hatch and spiral staircase vent to 
the outside was set to a 5.28 m’ rectangle. 


Volume from Level 1.7 to Dome 
Floor Openings Connect to 
Level 1.6 Staircases and Hatch 


Main 
Staircase 


Spiral 


Hallway Staircase 


Fire Room + 
Doorway 


Indicates Connection to Outside 
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Figure 7.5: Level 1.600 Figure 7.6: CFAST Model Block Diagram 
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The remainder of the 1.600 level included all of this level except for rooms 1.602 and 1.609 
which were closed off. This resulted in a compartment volume of 288 m°. The resulting 
compartment dimensions using a square floor were 7.92 m x 7.92 m x 4.60 m. All room 
surfaces were set to 0.50 m thick concrete. Vent connections were established between this 
compartment and the other three compartments on the level as well as between the main 
staircase and hatch. 


The portion of the HDR facility above the 1.600 level was modeled as a single compartment. 
The volume of this compartment was obtained by summing up the volume of all compartments 
above the 1.600 level excluding rooms 1.702 and 1.706. This yielded a total volume of 7215 m° 
and preserving the full HDR height resulted in a compartment with dimensions of 14.37 m x 
14.37 m x 39.95 m. All wall surfaces for this compartment were set to 0.5 m thick concrete 
which is not a true representation for the steel shell dome. The level 1.600 hatch and staircase 
openings to the were connected to this compartment. 


There was one final compartment that was incorporated into the model. This compartment was 
added to allow flow to return from the dome, down the spiral staircase, and then down to the 
1.400 level at room 1.402. Even though the spiral staircase hatch to the 1.500 level was closed, 
there did exist numerous floor penetration in the form of pipe channels, cable trays, and other 
conduits that connected the 1.600 level to the 1.500 level and the 1.400 level and these need to 
be included to properly model the HDR facility. The room connection list [1] was consulted and 
the area of available connections from the 1.400 level to the 1.600 level (excluding hatches and 
staircases) was totaled. This yielded a flow area of 9 m*. Since CFAST would not initialize 
successfully with a completely enclosed geometry, a vent to the outside was created in this 
compartment which was closed after the 100 s null transient portion of the model. The resultant 
compartment started at an elevation of 0.3 m and extended vertically for 13.2 m, 2.4 m above the 
floor of the 1.600 level. The 1.600 level compartment connected with a 3.0 m x 2.4 m doorway, 
room 1.402 connected with a 1.8 m x 3.2 m doorway, and the vent to the outside was 0.05 m 
wide that extended from the ceiling of room 1.402 to the floor of the 1.600 level. 


A block diagram of the wood crib model is shown in Figure 7.6 on page 7-7. The final CFAST 
model contained eleven compartments on three levels, the 1.400 level and part of the 1.500 
level, the 1.600 level, and levels above the 1.600 level. These compartments were connected by 
a total of eighteen vents of which five were vertical flow vents. There was one vent to the 
outside which was closed before the start of the fire. Shown on the next page below are the 
CFAST input cards for this model. 
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8 CFAST RESULTS AND COMPARISONS WITH HDR DATA 


This section compares the CFAST predictions using the model described in Section 7 with 
measured data from the actual tests for each of the three wood crib tests. Instrument descriptions 
and locations can be found in Section 3 of this report. The text and figures that follow will 
compare predictions made for the three wood crib tests with CFAST predictions to examine the 
effects of fire power and fuel availability on CFAST’s predictive capabilities. In each of the 
figures that follow, the HDR sensor is identified by its instrument number and its axial position 
in terms of the HDR coordinate. 


As this section discusses model comparisons with data, it is important to define what the authors 
consider a good versus a poor comparison. CFAST is designed to be a quickly executed 
engineering tool with a relatively small learning curve. Therefore, CFAST cannot be expected 
to make exact or near exact predictions to the data especially given the complexity of the HDR 
facility. However, it can be expected that CFAST should reproduce the same regional and 
transient trends as seen in the data, predict with reasonable accuracy the significant phenomena 
of the experiments, and not introduce significant non-existing phenomena. 


The first set of figures, Figures 8.1 through 8.3, shows the pyrolysis rates used as inputs to 
CFAST vs. the pyrolysis rates as measured by the digital scales. It is clear from these figures 
that the pyrolysis rate utilized by CFAST closely matches the data, both in terms of maximal 
values and net mass input into the fire room. For the three tests the deviation from net mass 
input is 3.5%, 1.0%, and 0.4% respectively. 
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Figure 8.1: T51.16 Pyrolysis Rate 
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Figure 8.2: T51.17 Pyrolysis Rate 
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Figure 8.3: T51.18 Pyrolysis Rate 
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The next three figures, Figures 8.4 through 8.6, show the upper layer temperature in the fire 
room. The peak temperature predictions for all three tests are high with the maximum 
increasing with the number of wood cribs. The CFAST predictions agree surprisingly well 
during the heatup phase and match the transient temperature gradients of the fire room upper 
layer extremely well. However, the predicted cooldown is much to fast and thus leads to 
substantial non-conservative underpredictions. In other words, for the conditions examined, 
CFAST appears to be transferring energy appropriately to the fire rooms surfaces during the 
heatup phase, but is not releasing that energy correctly during the cooldown phase. CFAST 
results for test T51.18 show a puzzling behavior of the temperature increasing 80 ° C above the 
cooldown temperatures well after the fire was over As the only difference between the cases 
was the fuel availability, this indicates a possible problem within CFAST’s algorithms.. 


The three figures following Figures 8.4 through 8.6 show the lower layer temperatures in the fire 
room, Figures 8.7 through 8.9. As with the upper layer, CFAST overpredicts the lower layer 
temperature with the prediction worsening as fuel availability increases. Unlike the upper layer, 
CFAST continues to overpredict the temperature during the cooldown period after the fire. 
Lastly, during tests T51.17 and T51.18 there are plateaus and secondary peaks in the measured 
data which are not predicted by CFAST. This last observation most likely results partially from 
artifacts caused by the approximated smoothly varying pyrolysis rate given to CFAST as input 
and mostly from the inability of a zone model code to calculate local temperature behaviors 
resulting from detailed geometric aspects of a compartment. It is quite obvious from the 
comparison between data and CFAST predictions for all three wood crib tests that the agreement 
both with respect to transient characteristics as well as quantitative values of the temperature in 
the lower layer leaves much to be desired. On the other hand, the predicted temperatures are 
seemingly conservative. 
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Figure 8.4: T51.16 Fire Room Upper Layer Temperature 
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Figure 8.5: T51.17 Fire Room Upper Layer Temperature 
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Figure 8.6: T51.18 Fire Room Upper Layer Temperature 
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Figure 8.7: T51.16 Fire Room Lower Layer Temperature 
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Figure 8.8: T51.17 Fire Room Lower Layer Temperature 
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Figure 8.9: T51.18 Fire Room Lower Layer Temperature 


The two figures the next page, Figures 8.10 and 8.11, display the CFAST predicted fire room 
layer height for test T51.16 and the measured temperatures from the thermocouple rake located in 
the fire room doorway. The thermocouple temperatures indicate that the layer height in the fire 
room lies somewhere between the -0.1 m elevation and the +0.5 m elevation. CFAST predicts the 
layer height to be just above the floor of the fire room near -0.7 m. This is a significant 
underprediction of the layer height. Although it does explain the elevated temperatures predicted 
in the previous lower layer figures, as the almost nonexistent layer will experience a larger 
temperature rise from radiative heat. It does not, however, explain the upper layer temperature 
results. A significantly larger upper layer should tend towards underprediction due to the larger 
mass present in the layer. Results for the other two tests are equivalent. 


The next three figures, Figures 8.12 through 8.15, show the upper layer velocity of the hot gasses 
exiting the fire room. Note that the measured velocities are not zero at the beginning of the tests. 
This is due to naturally existing circulation patterns with the HDR facility. For each of the three 
wood crib tests, CFAST underpredicts the velocity by a factor of six. Accounting for the 
overprediction in the upper layer size by close to a factor of two, CFAST is still underpredicting 
the velocity by a factor of three. Also CFAST starts to predict the velocities reasonably well for 
the first five minutes of the fire, but then the predicted velocity drops suddenly to near zero before 
recovering slightly. This drop and recovery may indicative a possible instability within the code’s 
predictive capabilities. Lastly, at the end of test T51.18, CFAST predicts a slow rise in velocity, 
followed by a sharp oscillation. Though this echoes the hot layer temperature results, since the 
pyrolysis rate has been nearly zero for thirty minutes, it is completely nonphysical. 
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Figure 8.10: T51.16 CFAST Fire Room Layer Height 
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Figure 8.11: T51.16 Measured Fire Room Layer Height 
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Figure 8.12: T51.16 Fire Room Upper Layer Velocity 
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Figure 8.13: T51.17 Fire Room Upper Layer Velocity 
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Figure 8.14: T51.18 Fire Room Upper Layer Velocity 


CFAST predicted oxygen concentrations are shown with the corresponding measured data for the 
fire room hot layer in Figures 8.15 through 8.17. During the fire for tests T51.17 and T51.18, 
CFAST predicts very well the O2 concentrations measured in the upper layer of the fire room. 
However, the CFAST prediction of the O concentration for test T51.16 fails to match the data. 
This failure most likely results from the poor agreement of the layer prediction. A higher layer 
position during the end of the fire would have the effect of lowering the predicted O, 
concentration and improving the prediction. At the end of test T51.18, a nonphysical result 
shows up in the prediction. With almost all the fuel consumed by combustion, the O2 
concentration drops unexpectedly to zero for the last 20 minutes of the cooldown period. 


CFAST Results and Comparisons With HDR Data 8-9 


T51 WOOD CRIB FIRE REPORT NUMAFIRE:03-98 


NO 
NO 


CFAST 
20 | 
oi CG4642 
S. 18 +1.0m 
c 16 


Oxygen Concentratio 
iON FP DWOOON F 


1090 O8104200730 (8400150 250 970 GSO 0 
Time (min) 


Figure 8.15: T51.16 Fire Room Upper Layer O, Concentration 
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Figure 8.16: T51.17 Fire Room Upper Layer O, Concentration 
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Figure 8.17: T51.18 Fire Room Upper Layer O, Concentration 


The next set of figures, Figures 8.18 to 8.20, shows the CO,, CO, and unburned hydrocarbon 
concentrations in the fire room upper layer for test T51.18. CFAST CO, predictions for T51.16 
and T51.17 are equivalent to the predictive capability shown for T51.18, and for test T51.17 and 
T51.17 CO and unburned hydrocarbons are not a significant parameter due to the well ventilated 
nature of the fires. The CO) prediction by CFAST in Figure 3.18 is a factor of two lower than the 
measured concentration. However, since the layer height is underpredicted by CFAST, this low 
CO, prediction is to be expected. CO predictions during the fire in Figure 3.19 are also 
underpredicted, but by a much larger factor. Therefore, the layer height error cannot explain all 
of the CO deviation. Since CO production in CFAST relies mainly on the user input for the CO 
to CO, production ratio during ventilated burning, a large portion of the CO discrepancy probably 
results from a poorly chosen input value for the CO card. As is the value is chosen by the 
preprocessor when a wooden object fire was created, it suggests that better values might be 
required for that object definition as well. The unburned hydrocarbon concentration during the 
fire is slightly overpredicted by CFAST. This overprediction probably results from the 
underpredicted velocities in the fire room doorway. This would result in CFAST predicting the 
fire to be more oxygen starved than it was during the actual experiment and thus boosting the 
unburned hydrocarbon production. 


Each of the combustion product predictions in the fire room show a large concentration increase 
at the end of the fire by. The CO and CO; concentrations increase by a factor of four from their 
minimum post-fire values. The unburned hydrocarbon concentration increases by a factor of 
almost one thousand, to nearly 35% from almost 0%. This cannot be attributed to the hot layer 
collapsing back to the ceiling after the fire as the layer reduces by a factor of three whereas the 
unburned hydrocarbon concentration increases by nearly three orders of magnitude. It would 
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appear that CFAST is creating a source of fuel from nothing that is combusted as an 
underventilated fire. This hypothesis would explain the exotic behavior of the O) concentration 
behavior displayed in Figure 8.17, the temperature anomaly in Figure 8.6, and the velocity 
oscillations in Figure 8.14. 
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Figure 8.18: T51.18 Fire Room Upper Layer CO, Concentration 
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Figure 8.19: T51.18 Fire Room Upper Layer CO Concentration 
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Figure 8.20: T51.18 Fire Room Upper Layer Unburned Hydrocarbon Concentration 


Figures 8.21 and 8.22 on the following page display the predicted vs. measured velocities for 
test T51.17 in the main staircase and spiral staircase maintenance hatches at the 1.700 level. 
Results for tests T51.16 and T51.17 have a similar predictive quality. The CFAST predictions 
for the main staircase are below the measured values by a factor of three. This is similar to the 
predictive capability seen in the fire room velocities after taking into account the error in the 
layer height. This discrepancy is consistent with other velocity predictions. Surprisingly, the 
velocity prediction in the spiral staircase predictions agrees nearly exactly with the data. The 
prediction is somewhat below the data; however, since the measured value was located in the 
center of the hatch and CFAST’s prediction is for the average flow, one would expect the 
CFAST velocity to be slightly underpredicted. 


The final two figures, Figures 8.23 and 8.24, show the CO, concentrations and the extinction 
coefficients respectively in the dome for test T51.18. Again CFAST predictive abilities were 
equivalent for the.other wood crib tests for these quantities. The CO, concentration in the dome 
is underpredicted by a factor of two by the end of the cooldown period. This is the opposite of 
what was expected. Since the HDR as modeled in CFAST contains less free air volume than the 
actual facility, one would expect CFAST to overpredict the CO, concentration no underpredict 
it. The extinction coefficient prediction during the cooldown period is quite accurate; however, 
the prediction during the heatup phase of the fire is not. CFAST fails to predict the sharp peak 
in the extinction coefficient associated with the highly underventilated portion and does not 
calculate any decrease in the optical density at the end of the fire. The latter results from 
CFAST not accounting for the aerosol behavior of smoke, i.e. settling, agglomeration, etc. 
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Figure 8.21: T51.17 Level 1.700 Main Staircase Maintenance Hatch Velocity 
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Figure 8.22: T51.17 Level 1.700 Spiral Staircase Maintenance Hatch Velocity 
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Figure 8.23: T51.18 Dome CO, Concentration 
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Figure 8.24: T51.18 Dome Extinction Coefficient 
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9 CFAST OBSERVATIONS AND COMMENTS 


This section documents the general usability impressions the authors had while using CFAST for 
simulating the HDR: T51 wood crib fire experiments. The first part addresses the 
accomplishments by and with CFAST over the spectrum of experiments covered during this 
validation effort. The second part covers differences and deviations in comparison with the data, 
behavior during code execution, performance of the implemented models, comments about the 
documentation, and comments about the capabilities and limitations of pre- and post-processors. 
The comments are made with respect to performance based trends as they impose requirements 
in terms of quality. The chapter closes with suggestions for continued validation. 


9.1 Accomplishments with CFAST 


The following accomplishments were achieved with CFAST during the validation efforts using 
the wood crib fire tests T51.16, T51.17, and T51.18: 


¢ CFAST v3.1.1 was capable of modeling the HDR facility as a completely enclosed volume 
after a brief prefire period with an opening to the outside. 

¢ CFAST predicted upper layer temperatures in the fire room agree very well during the 
heatup phase with the data for all three experiments. The maximum temperature predicted 
by CFAST tends to be conservative with increasing fire power. 

¢ Oxygen concentration in the fire room during the fire is well predicted. 

¢ CO, concentration in the fire room is reasonably well predicted once layer height errors are 
accounted for. 

¢ Unlike the predictions for the gas fire tests, CFAST computed O, concentration in the fire 
room returned to ambient levels after the fire ended. 

¢ During the underventilated portion of test T51.18, CFAST performed well in calculating the 
percentage of unburned hydrocarbons. 


9.2 User Observations 


This subsection documents various aspects of code usage, applications, and executions that were 
observed by the authors during this project. 


9.2.1 Observations During Code Execution 

These are some observations made regarding the version of CFAST used for this report: 

¢ The added runtime feature of alerting the code user when time steps become very small for 
extended periods is a useful one. 


e The new postprocessing utility, REPORTSS, which allows for the output of CFAST results 
in a spreadsheet compatible format is also very useful. 
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9.2.2 Physical Models 
The following are some observed deficiencies in the CFAST physical models: 


¢ The O2 card is not implemented properly within the code as is shown in Appendix B. 
Furthermore, the code documentation makes no mention that the O2 card has 
implementation difficulties even though this has been the case for many years. 

¢ The CFAST model for the wood crib fires calculated essentially no flow in or out of the fire 
room when in fact substantial flows existed. 

* CFAST predicts the lower layer height much too low, which leads to a substantial 
overprediction of the upper layer thickness. 

¢ CFAST does not appear to compute any settling of particulate matter which acts to reduce 
the optical density and improve visibility in the far field from a fire and in the long term 
cooldown after a fire. 

¢ CO, levels in the facility were predicted to be much lower than indicated by the data. 

¢ During the underventilated portion of test T51.18, CO production did not substantially 
increase as it should during an underventilated fire. 

¢ At the end of test T51.18, CFAST calculates a significant amount of combustion products, in 
particular unburned hydrocarbons, even though no more fuel is available in the fire room. 


9.3 CFAST Status 


Based on all of the experiences with the CFAST code at UMCP in comparison with the HDR 
data, code performance, utilization and applicability have advanced over the last years. In fact 
as documented by many contributions at national and international conferences, the CFAST code 
has left the realm of exploratory research and entered into the competitive, highly demanding 
industrial market and regulatory decision making. 


For many aspects the CFAST predictions satisfy the stringent requirements. For others, 
however, it dramatically fails in a rather erratic manner and one has to wonder how the favorable 
agreements with data may change once the major remaining discrepancies have been resolved. 


The current efforts of UMCP, Department of Materials and Nuclear Engineering, were not 
designed to solve these broader aspects of balanced code behavior. They rather contribute 
towards the identification of remaining issues. 


However, it seems appropriate to point out that additional efforts are needed to maintain 


CFAST’s current position and impact. This especially in view of the fast evolution of 
performance based methodologies. 
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9.4 Suggestions for Further Validation 


One fact that is clear from the work performed for this report is that further work is needed in 
resolving the exact nature of the O2 card failure. Further test of this card along with a detailed 
examination of its use within the source code may resolve this problem. 


Only one of the wood crib tests became completely underventilated. CFAST results for this test 
were mixed. While CFAST did predict well the O. and unburned hydrocarbon levels, it did not do 
as well predicting other combustion products or the optical density. Further validation efforts 
with underventilated fires are needed especially considering that it is the high smoke and CO 
production by such fires that pose some of the greatest risks to human life and intervention 
efforts. 


It is suggested that the CFAST code developer group consistently evaluate some of the observed 
deficiencies. The UMCP input models will be made available electronically. Lastly, it is clear that 
the default values for CO and soot/smoke production used in the predefined object descriptions 
do not function well in modeling real fires. More effort is needed in validating and developing 
reliable values for use in modeling real fires. 
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APPENDIX A: INPUT FILES 


This appendix gives the complete CFAST input files for all computations documented in this 
report. The filename nomenclature isT51XX-C.INP where XX indicates the test number in the 
wood crib test series. 


T5116-C.INP:CFAST input file for Test T51.16 (5 cribs) 


VERSN 3T51-16: 9 Room+Dome+Level 1.4 Connect 

TIMES 5500 100 50 0 0 

DUMPR T5116-C.HIS 

STPMAX 1.00000 

TAMB 288.000 101300. 0.000000 

EAMB 288.000 101300. 0.000000 

HI/F OnZo0nnO  COOmO. OOO 0. C0ORO COOMIa SROO wl OO CLL FOO mie: SOO RL 62 25.0 10!53'00 
Weir OS Oe 950mm. 400) 4 350 eel. 800) 7.920 42300 45330, 9.670 242370 3.000 
Die tee sO FON COOMA MO Ole. H20hnO + SO OMea? aoe OMe 2 a5 0 82 A720 8 UO S67 O04 23703. 0.00 
HetGh eae SUM eoomi OO dLe LOOMS 500) "4,600. 5.0505 S050) 52050 34.950 13.200 

CEILI FIRECEIL YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
WALLS YTONG250 YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
FLOOR FIRE_FLR CONCR100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 


BVENGS 2 2P Ae LROLOR W975" 10000 

CVENE, £2 il O00) 2.0005 2.000 L000, “22000 S000 “2000 “15000 A000 1.000: T2000 
VENT, 253) 10 15.800) 2.485) 0)..000 

CVENTS 2a see 0007" LiIOO0> 1000) “22000 “Hs000) O00 S000. 000 sZ000 L000) a.000 
HVENE 39 49° 1) 47.300) 07.500 01.000 

CVENT oe Cee 5.000) E000. T6000) 2000) T55000) SEI000 HE000" VET000! 2000) 1.000 “13000 
HVENT, 3) 5) 2) 2.000) 0.500) (01.000 

CVENTA SS) sone 000) U7 000)5 L000) 1000) YATO00r =25000' “2000 x0do: 1.000 000 (15000 
HVENE 465) 1) 22700) 3.2007 1.200 

CVENT, 7405 >a) 1) 17000" 1000} £000) 2.000: “1000 «2000 W000 “Lr000) 1000' 2.000 -1.:000 
HVENT> 5) Li i) 0 '800)° 3'.5007  0).300 

CVENI >) dale ty 1000) 000% 12.000) ~1000F 22000) T5000 “O00 25000 22000 “12000 120100 
HVE » 6s 70 Ty 3S000 4.250 0/000 

CVENi Orie sOOCn a,.000 St OOOR tT O00) 12000) 12000 1,000 L000 12000 1.000 1.2000 
HAVENT 6 8 2 3.000 4.250 0.000 

CvENi Son on aee O00 dt O00 OOO T0000 15000 12000 2.000 1.000 2:5000 * 1.000 —-1-..000 
BVENE (67 9) 132000) 4.250 (0.000 

CVENin | Comoe OOO) SL A000 F000 000 2000" SasC00 1.000 1.000 “L000 2.000) F000 
HVENT 7 8 1 3.000 4.250 0.000 

CVENE ee oeeetes b OOOR IT O00 a 000) 12000" 22000) 12000 2.000 1.000 17.000" 1.000) 4.000 
BVENG ete O030 9:.800 3.3100 

CVEN Tae ee O00 1 O-000 0000 0-000 O2000 F000 0-000 0.000 0.000" 02000 0-000 
PVENG Guise 3. 000° 2.400 *0).:000 

CVENCe Gpiieee = L000 i 0008 O00: “L000 F.000 3000 10001 7000) S3000F” VI000)S 1.2000 
VVENT 7 3 4.54000 2 

VVENT 8 4 5.75000 2 

VVENT 10 7 4.54000 2 

VVENE 103 1896.97000 2 

VVENF 10° 97) 5. 28000" 2 

CHEMI 13.800 20.000 0.000 2.721E+007 293.000 493.000 0.200 

LFBO 1 

LFBT 2 

CJET OFF 


FPOS Aso d.825.9 0.000 
FTIME 100.00 700.00 1000.00 1126.00 1306.00 1480.00 1720.00 2260.00 3100.00 4900.00 


FHIGH 0.00000 0.00000 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.00000 
FAREA 0.00000 0.00000 0.70900 0.70900 0.70900 0.70900 0.70900 0.70900 0.70900 0.70900 0.00000 
FMASS 0.00000 0.00000 0.00926 0.04260 0.03705 0.01806 0.02454 0.01158 0.00232 0.00116 0.00000 
FQDOT 0.00000 0.00000 252000. 1159200. 1008000. 491400. 667800. 315000. 63000.0 31500.0 0.00000 
HCR 0.00000 0.00000 0.21050 0.21050 0.21050 0.21050 0.21050 O° 21050 0.21050 0.21050 0.00000 
oD 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
co 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
#GRAPHICS ON 

DEVICE 1 

WINDOW Or he OF 1279. 1023". 4095. 
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T5117-C.INP:CFAST input file for Test T51.17 (7 cribs) 


VERSN 3T51-17: 9 Room+Dome+Level 1.4 Connect 

TIMES 5500 100 50 0 0 

DUMPR T5117-C.HIS 

STPMAX 1.00000 

TAMB 288.000 101300. 0.000000 

EAMB 288.000 101300. 0.000000 

HI/F 0.250: 0.000) 10,7000 G.O00M 02000 RI LOOMS LOOM OO ees O OM Gn 5 OneeOc a0) 

WIDTH 2:..950 4.45950 54.300 4.3830 .17800 972920\9 4,300. 4° 330.5 9.6 7,0N142370) 30000: 

DEPTH 3).650! 27800527750) 22720) 76.350) 72920) 9'257500 2/20 o 670) 4 23 Owe 0 0.0 

HETGH 2.750 2.485 12.2100) 11.100, 32500, 4.600., 5,050" 52050) 52050349502 133200 

CEILI FIRECEIL YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
WALLS YTONG250 YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
FLOOR FIRE_FLR CONCR100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 


HVENT (15°29 Via OL 0) S975 On 000 

CVENT, 1 2 1) 100005270005 L000) 13000 ~1 C00 120007 14 3000 MeLE COO, HA RO0O a 0008 eE000 
HVENT® 62 35) eel S00 24859 0000 

CVENT 2 3 de als000) 1.0003) 12000) 12000 12000" 2000) 42000 2a C0C Ty EO00 000m e000 
HVENT 3 4 1 473007 01500) 0000 

CVENT 3 4. 1 2.000 LA0009 2-000) “2000 1000) LACO0 L000 Fes 000) Ses O 00s aOC0 me 000 
HVENT 3 5) - 1) 220009025008 30/000 

CVENT 3 5 d 1.000 (§1.:000,.1..000) 4.3000, 1000) 51 7000/5 1 S000 R000) Be es00 0.6 L000 aango 0 
HVENT 4° 5 12700 RRS sL00) & 1200 

CVENT 4 5 1 2.000 21.000 TROOO MRIS 000M 12000 MAO OOr 0008 T2000 e000) 4000 1000 
HVENT 5 id» 1.128005 3.500! 05800 

CVENT 5 12> dy 012000) © 2.0005) T4000) (a OOO A OO OP O00) LAC OORT A FONG 000 RL 000m 000 
HVENT ‘6°47 925535, 000) -4.:250'5) 05 000 

CVENT 67 Jailed O00r L000 Ral, O00) 5 L000 sO00F SOO 0s 552 .0:0,0 Ril C.00 ies 6) 0/0) eaten 00s aa 0.0)0) 
HVENT’ 6 87 2933;10007) 435250550)..000 

CVENT 6 8 2 2.000 20005511000, 1000" “S000 GaO00) Ss 000 | F000, 000g 000 Sks000 
HAVENT 6 -9 dL? 32,000) 745250! 70.000 

CVENT 6, 9 1 2.000, 22000 2,000) 1.000) GEx000 LERO00 M000) TsO 007 000)  eis000) e000 
HVENT 7 8° ASSO 08 V4E 250 0000 

CVENT 7 8 1 1.000) 12000 1.000 1.000 HROO0) Sie000 M000" E000) W000) Ti000) Maz0o0 
AVENE 22 07030) 92800 wes s00 

CVENT 11 12 1 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
HVENT 6 117 -1 3.000 2.400 0.000 

CVENT’ ‘6 Ji 21 2/000) -22000,75 P2000) LI000, VaesO00F HO O0: Ta O00> TIRC00" MAO00, SOOO” Simo 
VVENT 7 3 4.54000 2 

VVENT 8 4 5.75000 2 

VVENT 10 7 4.54000 2 

VVENT 10 8 6.97000 2 

VVENT 10) 99 5.280002 

CHEMI 13.800 20.000 0.000 2.721E+007 293.000 493.000 0.200 

LFBO 1 

LFBT 2 

CJET OFF 


FPOS 12475 17825: 7105000 
FTIME 100.00 700.00 1000.00 1144.00 1300.00 1780.00 2380.00 ST00C00 27 5500-00 


FHIGH 0.00000 0.00000 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.00000 
FAREA 0.00000 0.00000 1.10300 1.10300 1.10300 1.10300 1.10300 1.10300 1.10300 0.00000 
FMASS _ 0.00000 0.00000 0.01389 0.05326 0.05372 0.04168 0.01389 0.00463 0.00139 0.00000 
FQDOT 0.00000 0.00000 378000. 1449000. 1461600. 1134000. 378000. 126000. 37800.0 0.00000 
HCR 0.00000 0.00000 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.00000 
OD 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
ee) 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
#GRAPHICS ON 

DEVICE 1 

WINDOW QO. Of On L279) 10235 40955 


Appendix A: Input Files A-2 
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T5118-C.INP:CFAST input file for Test T51.18 (11 cribs) 


VERSN 3T51-18: 9 Room+Dome+Level 1.4 Connect 
TIMES 5500 100 50 0 0 

DUMPR T5118-C.HIS 

STPMAX 1.00000 . 

TAMB 288.000 101300. 0.000000 

EAMB 288.000 101300. 0.000000 


HI/F OeZ50 NF OF0008 OF000) 0.0005) OF 000 11100 PL rO0: 117100) TL 100) 26.150) 10:.3.00 
WIDTH sn 2;. JoOmsaraI oO 4523.00 4° 330' 800) 7.920 4.300 4.330 9.670 14.370 3.000 
DEEL OsUmeL OOO me. 7/50) 82720 96.350 997 3920) 22750" 2.720 §9.670 145370 632000 
Hoth eco SU econ LOOM dT 100 93.5005 74.600) —5-050' 5.050) 5.050 343950: 13.200 


CEILI FIRECEIL YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
WALLS YTONG250 YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
FLOOR FIRE_FLR CONCR100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCR050 


VENTS Se 2ee eee OLO EST 5 0.000 

CVEN Ri ea eee OOO b OOO R000 W000 12000) L000 1.0007 1.000 L000 41.000 21.000 
HVEND 2 Sede) L800 2-485 7 0)..0.00 

CVENESS2eeS a. L000 15000 L000 L20009 LACO0” 250008 E2000) TA000) LICOOr 13000)" 12000 
HVENE 3) 4250s) 453007 02500) 02000 

CVEN Ie ema OOO OOO 1 1000 1 000% 1 S000 130008 L000. 17000) LIC00F L000) 17000 
HVENT 3) (55a 20002 025007 01.000 

CVENT se See O00 met O00 1 OOO TAO O0) 000M? OOO REL OOO E0007 000 L000K ak. 000 
HVENT 4° Sl 25700) (3). 100M LOO 

CVENT Asa OOO de, OOO T0001 S000" FL000F 2000S LI000F L000, e000)” T2000)” 2.000 
HVEND) 5) Leet 12800 63:.500  10).300 

CVEN TS eee OOO O00 Reds OOO Is. O00 4.000 il 0002000) L000 000 Mal 000) 1-000 
HVENT “6 Feeds 35000 4.250 <0:.000 

CVENE Om ee OO Onl O00 Ne  OOORmHy OO 0R a. 000 teas 0007 000 St OC OM. 0007-000 TEis000 
HVENT 565) 89945 35000) 4.250) 0). 000 

CVEN ET SeOmG mere O00) Eb OOO. 000 MEE O00 1k 000M 000 Pe. 000 L000 1: 000TeL000 2.000 
HVENT | 6 9 S5ee 3'2000) 4.250950. 000 

CVEN TO een 2) 2000) O00 Mas 000 O00 e000) 00077 1.000 > 1.0007 Ay.000> 5. 000 7 2. 000 
HVENE, SSeS 2000) 4.250 180. O00 

CVENT omen O00 me OO0R. OO0ME. O00. 000. 000". 000" 1.000 —1)..000 921.1000 1.000 
HVENT, Li tere 05030) 9800773). 3:00 

CVENT 11 12 1 1.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
HVENE 6) 2290i 33000 2... 400P80).000 

CVEN Ome en O00 sed O00 tel 0 OOM. OOO. O00 a. 000 1: 000 1-000 1.000, 1.0007 1.000 
VVENT 7 3 4.54000 2 

VVENT 8 4 5.75000 2 

VVENT 10 7 4.54000 2 

VVENT 10 8 6.97000 2 

VVENT 10 9 5.28000 2 

CHEMI 13.800 0.000 0.000 2.721E+007 293.000 493.000 0.200 

LFBO 1 

LFBT 2 

CJET OFF 


FPOS TEATS. 8255 70.000 
FTIME 220.00 700.00 1060.00 1300.00 1600.00 1900.00 2500.00 3700.00 4900.00 5500.00 


FHIGH 0.00000 0.00000 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.00000 
FAREA 0.00000 0.00000 1.84300 1.84300 1.84300 1.84300 1.84300 1.84300 1.84300 1.84300 0.00000 
FMASS 0.00000 0.00000 OROT389 0.08567 OL 06715 0.04862 0.00926 0.00463 0.00232 0.00116 0.00000 
FQDOT 0.00000 0.00000 378000. 2331000. 1827000. 1323000. 252000. 126000. 63000.0 31500)-'0 0.00000 
HCR 0.00000 0.00000 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.00000 
OD 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
SO 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
#GRAPHICS ON 

DEVICE 1 

WINDOW O. 0. QO. 22793 1023\. 4095. 


Appendix A: Input Files A-3 
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APPENDIX B: CFAST O, CARD 


This appendix discusses the difficulties encountered when using the CFAST O2 card. The input 
file for T51.16 was modified to add the O2 card. These modifications consisted of respecfying 
the heat of combustion and pyrolysis rate back to their unmodified values. The modified T51.16 
input file is shown below: 


VERSN 3T51-16: 9 Room+Dome+Level 1.4 Connect 

TIMES 5500 100 50 0 0 

DUMPR T5116-D.HIS 

STPMAX 1.00000 

TAMB 288.000 101300. 0.000000 

EAMB 288.000 101300. 0.000000 

HI/F 0-250) OOOO 000M O2000 0 SO00MEs VOOM SOOM LOOT LOO Le. E50. 1053.00 

WEP 950 9 4,950" 43008 423305 1.3800 9772920" 4.300" 4.5330) 92670 14.370" 35000 

DEP LHe ODO aa sOOF een oO) 24/720 6.550 75920 8825750 92.720) 97670 145370) 3.000 

HeiGHeea (50n 2e485 ii 100) Lielooe 35500 647600 ~52050 52050) S2050) 342950) 13,.42100 

CEILI FIRECEIL YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
WALLS YTONG250 YTONG100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 
FLOOR FIRE_FLR CONCR100 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 CONCRO50 


VEN Dee 2s 010 1975) 705.000 

CVENT 1 2 1 COOREE O00 Tet 000) e000 2000 000) T0007 2000) 10005 1.0008) 1.000 
HVENG 2) 3 e800" 25485, 70.5000 

CVENT Eee St. O00 ae 000" 000 L000. £000 1 000F L000" L000 000 1000" as. 000 
HVENE 3 4 1 4.300 0.5007) 0.000 

CyENE WS 4 000000 1000 12 000 1.000 12 000F 1000" POCO 1000 “ak 0007 “2-000 
PVENT OS Lo 2 000 T0500 (0-000 

CVEN Teo ee O00 OOO 000 L000) L000 SPreCOr M000" L000 = 12000 1.0005 21-000 
VEN Ge eS Me a LOOe eS LOO 1100 

CVE Tee LOO Oma OOO L000 at S000 000 27000 ~ 40007 1.000 2.000 1.000" 1.000 
BVENT, 5S ik 2 1.800 7 3.5500 »0).300 

CVENT 55 tt O00) 2000 S000) V1 5000 = C00 =LROC0F L000 L000) M22 000" TF000 7 1.000 
HVENL. 6 7-2) 935- 000) 74.250 =0'000 

CVEND Ones bn OO0 O00 11000 1000) 1000, 12000; 12000 120007 21.000 13000 1.000 
EVENT Gass. i 3.0005 4.250 0.000 

CVENTenGO eS te. OCOM 12000" 1.000 25000) 1000 120007 1.000" 1.000 22000 17000) 2.000 
HVENE, 46-9 15 35000492508 05000 

GVENEF 1659 97 Se (000 FeiIeO00R 1 1000 ~ 1000.7 12000" “15000 —12000 120007 W000 — 150007 1.000 
HVENT 7 8. a 3 000) 42505 10...000 

COVEN aS, Leh. O00 is, OOOe E000 aEt000 23000" 1.000 71.000) 11.000 712000" L2000, “£2000 
AVENE Tis 2 1 On030, S9-.3005 3.300 

CVENP 1 12) i L000F W000] 02000" (050008 02000) “O2000"* 05000 §O5000' 00001" (02000 * 0) 000 
HVENT 611 1 3.000 2.400 0.000 

CVENIES Gel SA ICOOM LOCOS E000 000 MI000r E000. LEO00F LAO00 SP l000r 0007 17000 
VVENT 7 3 4.54000 2 

VVENT ~§8' 4 ~5.75000° @2 

VVENT 10 7 4.54000 2 

VVENT 10 8 6.97000 2 

VWVENT 10 9 5.28000 2 

CHEMI 29.800 20.000 0.000 1.262E+007 293.000 493.000 0.200 

LFBO 1 

LFBT 2 

CJET OFF 


FPOS IS475 5) 128259 012000 
FTIME 100.00 700.00 1000.00 1126.00 1306.00 1480.00 1720.00 2260.00 3100.00 4900.00 


FHIGH 0.00000 0.00000 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.28100 0.00000 
FAREA 0.00000 0.00000 0.70900 0.70900 0.70900 0.70900 0.70900 0.70900 0.70900 0.70900 0.00000 
FMASS 0.00000 0.00000 0.02000 0.09200 0.08000 0.03900 0.05300 0.02500 0.00500 0.00250 0.00000 
FQDOT 0.00000 0.00000 252000. 1159200. 1008000. 491400. 667800. 315000. 63000.0 31500.0 0.00000 
HCR 0.00000 0.00000 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.21050 0.00000 
02 0.00000 0.00000 1.40400 1.40400 1.40400 1.40400 1.40400 1.40400 1.40400 1.40400 0.00000 
OD 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
co 0.00000 0.00000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 0.03000 
#GRAPHICS ON 

DEVICE 1 

WINDOW QO. 0. 0.9 2279'11.023.e5 4095" 


The three figures below, Figures B.1 through B.3, show the calculated fire power, upper layer 
temperatures in the fire room, and the upper layer velocity from the fire room to the hallway. 


Appendix B: CFAST O2 Card B-1 
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Figure B.2: Fire Room Upper Layer Temperature 


Appendix B: CFAST O2 Card B-2 
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CFAST 


CV4640 
+0.96 m 


Velocity (m/s) 


9 bu 
10680 2910 620 93044405550 (26095700 '80 0.90 
Time (min) 


Figure B.3: Fire Room Upper Layer Velocity 


Figure B.1 and B.3 show that with the O2 card activated the fire generates the same power and 
mass flow rate as with the O2 card deactivated. Since fire power and exit mass flow from the 
fire room determines the heat input into the fire room, if they are equivalent between the two 
cases, one would expect the fire room temperatures to also be equivalent. Figure B.2 shows that 
this is not the case. With the O2 card turned on, the upper layer temperature in the fire room is 
over 200 °C higher than with the O2 card turned off. 


The next two figures, Figures B.4 and B.5, show the upper layer O) and CO, concentrations in 
the fire room. Since the empirical formula of wood is Co9sH240, each mole of wood empirical 
formula requires an additional 0.55 moles of O2 for complete combustion. Therefore, one would 
expect a decrease in the fire room O, concentration to be calculated by CFAST. CFAST, 
however, calculates a slight increase in the O2 concentration. Further more, since the same 
amount of carbon is being combusted, one would expect to see similar CO, concentrations 
between the O2 card and the no O2 card predictions. With the O2 card activated, CFAST 
calculates almost no CO; production. 


It is clear from these few comparisons that with the O2 card enabled difficulties arise due to its 
implementation and its use should be avoided if CFAST is to calculate reasonable results. 
Based one these observations it was decided to remove the oxygen from the wood when 
specifying the fuel in CFAST. 


Appendix B: CFAST O2 Card B-3 


T51 WOOD CRIB FIRE REPORT NUMAFIRE:03-98 


Oxygen Concentration (v/o) 


Carbon Dioxide Concentration (v/o) 


10480 10° 20%300,4095016070 680 70 
Time (min) 


Figure B.4: Fire Room Upper Layer O, Concentration 


CFAST 


CG4643 
+1.0m 


10 10 20 30 40 50 60 70 80 90 
Time (min) 


rO ja NWR AHDN OW 


Figure B.5: Fire Room Upper Layer CO, Concentration 


Appendix B: CFAST O2 Card B-4 
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